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SECTION  I 


INTRODUCTION 

The  problem  of  water  ingestion  into  engines  arises  both  when  water 
is  injected  deliberately  to  increase  the  thrust  output  and  when  water 
enters  an  engine  accidentally,  for  example  during  rain  or  operation 
from  a rough  runway  with  puddles  of  water  on  it.  A complete  analysis 
of  the  problem  of  water  ingestion  requires  consideration  of  an  installed 
engine  including  engine  location  with  respect  to  the  airframe  and  also 
the  wheels.  Thus,  one  has  to  establish  changes  in  the  performance  and 
stability  of  the  following. 

a.  inlet 

b.  fan  and  fan  duct 

c.  compressor 

d.  diffuser 

e.  combustor  and 

f.  afterburner. 

In  addition,  one  has  to  take  the  effects  of  water  ingestion  into  account 
in  engine-airframe  integration.  Among  the  components,  while  each  of 
those  referred  to  above  is  of  importance,  it  is  clear  that  the  compressor 
performance  and  stability  are  crucial  to  the  establishment  of  overall 
performance.  The  investigation  reported  on  here  is  an  extention  of  an 
earlier  study  (Ref.  1)  on  axial  flow  compressor  operation  with  air-water 
mixture  flow. 
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I 

1.1.  Outline  of  the  Report 

The  three  principal  topics  discussed  in  the  report  are  (1)  analysis 
of  selected  compressors  with  water  ingestion  and  evaporation,  (2)  engine 
instability  simulation  with  water  ingestion,  (3)  further  development  of 

the  axial  flow  compressor  performance  calculation  program  for  two  < 

1 

phase  flow.  The  later  includes  a discussion  of  the  method  of  incor- 
porating variable  composition  and  properties  of  the  working  medium  into 

the  performance  calculation  procedure.  The  three  topics  are  discussed  1 

in  the  following  three  sections  of  the  Report. 

j 

The  principal  objective  of  the  investigation  was  the  establishment  i 

of  the  major  areas  in  which  further  research  is  required  in  axial  flow  , 

compressors  when  operating  with  air-water  mixture  flow.  The  preliminary  , 

studies  reported  on  here  are  intended  to  identify  such  outstanding  pro- 
blems and  a short  discusion  on  them  is  provided  in  the  last  chapter 
entitled  "Discussion",  along  with  conclusions  from  investigations  on 
two  selected  compressors. 


SECTION  II 


1 

ANALYSIS  OF  SELECTED  COMPRESSORS 

When  water  is  ingested  into  a compressor,  two  general  questions 
can  be  posed  as  follows. 

(i)  Suppose  the  water  evaporates  at  a certain  station  in  the  com- 
pressor. What  is  the  effect  on  the  remaining  stages  of  the  compressor? 

( i i ) Suppose  again  that  the  water  evaporates  at  a certain  station  ■ 

in  the  compressor.  If  we  desire  to  bring  the  performance  of  the  stage 

following  that  station  to  the  values  that  existed  when  operating  with  ' 

air  alone,  all  other  conditions  being  equal,  is  it  possible  to  achieve 
this  with  bleeding  and/or  injection? 

Investigations  on  these  questions  follow. 

2.1  Influence  of  Blade  Loading 

When  water  evaporates  at  a certain  station  in  the  compressor,  it 
leads  to  the  following,  (i)  reduction  in  the  temperature  of  air;  (ii) 
change  in  the  composition  and  properties  of  the  flowing  medium;  (iii) 
change  in  the  corrected  speed  and  mass  flow  of  the  compressor  (iv)  change 
in  the  operating  points  of  the  stages  following  the  position  where  the 
water  evaporates. 

Consider  for  example  a 6-stage  compressor.  Suppose  we  assume  that 
the  initial  conditions  of  the  air-water  mixture  at  entry  to  the  compres-  i 

sor  are  such  that  under  a given  set  of  conditions,  the  water  undergoes  | 

evaporation  between  the  third  and  the  fourth  stages.  The  performance  of 
the  compressor  becomes  affected  in  the  first  three  stages  due  to  the 
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presence  of  water  ingested  with  air  and  because  of  the  upstream  effects 
induced  by  the  changes  in  the  successive  stages.  It  is  clear  that  the 
performance  of  stages  4,  5,  and  6 will  become  affected  because  of  the 
reduction  in  temperature  and  addition  of  steam  between  the  third  and 
fourth  stages.  The  performance  change  will  be  in  regard  to  the 
corrected  speed  and  mass  flow  and  pressure  ratio,  the  three  quantities 
that  make  up  an  operating  point.  Thus  the  overall  performance  of  the 
compressor  will  undergo  a change. 

In  presenting  the  modified  performance  of  the  compressor  for  a 
given  set  of  operating  conditions,  it  is  necessary  to  adopt  a uniform 
precedure.  The  one  adopted  here  is  the  standard  procedure  of  presenting 
the  compressor  performance  map  in  terms  of  the  mass  flow  parameter  at 
entry  to  the  compressor. 

2.1.1.  Problem  formulation 

The  following  assumptions  have  been  introduced. 

(1)  Water  droplets,  which  occupy  a negligible  volume  in  the 
flowing  medium  because  of  the  low  concentration  do  not  affect  the  com- 
pressor performance.  In  other  words,  a compressor  operating  with  air- 
water  mixture  operates  just  as  though  it  was  operating  with  air  alone. 

(2)  When  the  heat  content  of  air  and  water  are  right,  water  under- 
goes evaporation  instantaneously.  Consequently,  there  is  an  instantan- 
eous change  in  the  temperature  (reduction)  and  composition  (with  addition 
of  steam)  of  the  flowing  medium.  The  corrected  mass  flow  and  speed 
become  altered  and  therefore  the  output  of  the  stages  following  the 
evaporation  station  will  undergo  change. 
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(3)  The  mass  of  water  evaporated  is  assigned  a certain  value  in 
each  calculation. 

(4)  In  calculating  the  performance  of  the  compressor  after 
evaporation  of  water,  the  performance  of  various  stages  of  the  compressor 
is  based  on  ‘he  original  air  flow  data  (Refs.  2 and  3).  In  other  words, 
in  estimating  the  performance  of  a stage  with  air-steam  mixture,  we 
utilize  the  characteristic  of  the  stage  as  determined  with  air  as  the 
basis  while  accounting  for  changes  in  molecular  weight  and  ratio  of 
specific  heats. 

Utilizing  the  foregoing  assumptions,  the  performance  of  the  compres- 
sor with  evaporation  of  water  is  calculated  based  on  the  so-called  one- 
dimensional analysis.  In  that  analysis  one  assigns  flow  conditions  to  a 
certain  section  along  the  height  of  a blade  and  utilizes  performance 
characteristics  corresponding  to  that  section  of  the  compressor  to 
obtain  outlet  conditions  from  the  stage  at  that  section. 

It  should  be  pointed  out  here  that  the  following  investigation 
pertains  entirely  to  the  effects  of  evaporation  of  water  at  a chosen 
stage  in  the  compressor,  leading  to  a reduction  in  temperature  of  the 
fluid  locally,  the  exchange  of  mass  between  the  liquid  and  the  vapor 
phase  and  the  consequent  change  in  the  molecular  weight  and  ratio 
specific  heats  of  the  fluid  in  the  compressor.  Thus,  the  effects  of 
the  presence  of  water  in  some  of  the  cases  in  certain  stages  of  the 
compressor  and  the  induced  effects  of  changes  in  certain  stages  on 
other  upstream  stages  have  not  been  investigated  here.  Insofar  as  the 
presence  of  water  is  concerned,  in  view  of  the  small  volume  of  water 
ingested,  in  this  preliminary  analysis,  assumption  (i)  is  introduced 
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so  that  one  has  a reasonable  basis  for  carrying  out  the  calculation  of 
performance  when  there  is  air-water  vapor  mixture  in  the  compressor. 

It  may  also  be  pointed  out  that  the  first  attempt  at  modeling  air- 
water  droplet  mixture  flow  in  an  axial  compressor  is  presented  in  Ref. 

1 and  some  modifications  to  it  (to  account  for  different  droplet  sizes) 
is  presented  in  Section  4 of  this  Report.  No  calculation  schemes  based 
on  such  an  analysis  have  yet  been  evolved  and  the  models  themselves 
need  further  improvement. 

Regarding  upstream  effects  due  to  changes  in  a given  stage  of  the 
compressor,  no  attempt  at  evolving  a theory  has  been  undertaken  in  this 
Report. 

2.1.2.  Solution  procedure 

The  solution  procedure  and  a computer  program  (Purdue  University 
Compressor  Performance  Program,  PUCPP-001 ) are  described  in  Appendix  1 
and  Appendix  2 to  this  report. 

2.1.3.  Results  for  typical  compressors 

Calculations  of  changes  in  compressor  performance  have  been  preformed 
on  two  compressors,  namely 

(1)  N.6.T.E.  Reference  #109  Compressor  (Refs.  1 and  2)  and 

(2)  An  Allison  Research  Compressor  (6  stages)  for  which  data  were 
supplied  by  Detroit  Diesel  Division,  Indianapolis. 

In  both  cases,  we  have  investigated  the  effect  of  evaporation  of 
2.5%  by  weight  of  water  at  chosen  locations  in  the  two  compressors. 

There  is  no  special  significance  to  the  value  of  2.5%  by  weight  for  water 
that  is  supposed  to  undergo  evaporation. 
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It  may  be  pointed  out  here  that  the  N.G.T.E.  compressor  (Ref.  2) 
and  the  Allison  Research  Compressor  (Ref.  3)  differ  appreciably  in 
their  stage  characteristics  and  stage  loading.  Calculations  performed 
on  these  are  expected  to  show  the  effects  of  the  differences  in  those 
parameters. 

in  the  calculations  performed  on  the  NGTE  compressor,  the  entry 
conditions  have  been  assumed  to  be  as  follows. 

T = 485. OK 
°1 

PQ  = 1 atm.  (standard) 

The  temperature  of  485. OK  has  been  selected  for  air  as  follows:  suppose 
water,  at  a temperature  of  288. OK,  is  introduced  into  the  air  stream 
at  a mass  rate  of  5.0%  of  that  of  air;  then  at  the  air  stream  temper- 
ature of  485. OK,  all  of  the  water  would  be  able  to  undergo  evaporation. 
Thus,  if  the  air  temperature  is  485. OK  at  entry  to  the  first  stage  of 
the  compressor,  one  could  expect  complete  evaporation  of  2.5%  by 
weight  of  water,  which  is  the  case  investigated  here. 

Suppose  now  that  it  is  desired  to  obtain  the  evaporation  of  the 
same  amount  of  water  at  entry  to  the  second  stage.  Then  a larger  amount 
of  water  needs  to  be  ingested  at  entry  into  the  first  stage  such  that 
the  temperature  at  entry  to  the  second  stage  (equal  to  485. OK  plus  the 
rise  in  temperature  in  the  first  stage)  is  sufficient  to  evaporate  the 
required  amount  of  water. 

It  may  be  recalled  that  in  these  calculations  the  presence  of  water 
(prior  to  evaporation)  is  assumed  not  to  affect  the  performance  of  the 
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compressor.  Thus,  if  2.5%  by  weight  of  water  is  assumed  to  evaporate 
at  entry  to  the  4th  stage,  while  the  inlet  air  temperature  is  485. OK, 
even  though  this  would  require  the  presence  of  water  in  the  first 
three  stages  and  also  in  the  4th  stage  (in  the  latter  because  of  the 
presence  of  unevaporated  water),  the  performance  of  those  stages  is 
assumed  to  be  unaffected  by  the  presence  of  water. 

Calculations  have  been  performed  on  the  N.G.T.E.  compressor  at 
different  speeds  of  rotation  for  six  different  cases  corresponding  to 
the  same  mass  of  water,  namely  2.5%  by  weight,  undergoing  evaporation 
at  entry  to  each  of  the  six  stages  of  the  compressor  while  the  inlet 
temperature  is  held  constant  at  485. OK. 

It  is  important  to  point  out  that  a part  of  the  calculation  has 
been  based  on  extrapolation  of  the  stage  characteristics. 

In  the  case  of  the  Allison  compressor,  the  entry  conditions  have 
been  assumed  to  be  as  follows. 


T = 533. 2R 

Oi 


= 1 atm. 


(standard) 


The  temperature  of  533. 2R  corresponds  to  a chosen  test  condition  and 
therefore  an  arbitrarily  selected  value  for  the  purposes  of  the  cal- 
culation illustrated  here.  With  that  value  of  inlet  temperature,  it  is 
clear  that  by  adjusting  the  amount  of  water  that  is  ingested  into  the 
compressor  at  inlet  one  can  obtain  evaporation  of  2.5%  by  weight  of  water 
at  least  at  entry  to  certain  stages.  It  can  be  shown  easily,  considering 
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the  stage  temperature  rises  in  the  Allison  Research  Compressor,  that 
evaporation  of  2.5%  by  weight  of  water  can  occur  only  at  entrance  to 
the  4th,  5th  and  6th  stages  when  the  inlet  temperature  is  533. 2R. 

2. 1.3.1  The  results  of  calculations  on  the  N.G.T.E.  compressor  are 
presented  in  Figs.  2.1  - 2.3  for  three  speeds  of  rotation,  namely 
9,500,  8,000  and  6,000  R.P.M.,  respectively. 

In  each  of  the  Figures,  the  corrected  mass  flow  is  defined  as 
follows. 


corrected  mass  flow 


_ nv'o 


where  m is  in  lbm/sec,  0 = Tol/Tref>  <5  = p0i/pref»  Tref  = 288K»  and 
Pref  = 1 atm.  The  corrected  speed  is  defined  as  follows. 

N 

corrected  speed  = — 

where  N is  the  speed  of  rotation  per  minute. 

In  each  of  the  figures  the  performance  of  the  N.G.T.E.  compressor 
is  presented  in  the  following  manner 

(a)  Overall  performance  of  the  compressor  at  the  given  speed  for 
operation  with  air  only  at  an  inlet  temperature  of  288°K. 


7 _ , 


9 


(b)  Overall  performance  of  the  compressor  at  the  given  speed  for 
operation  with  air  only  at  an  inlet  temperature  of  485.00°K. 

(c)  Overall  performance  of  the  compressor  at  the  given  speed  when 


the  compressor  is  operated  with  inlet  air  temperature  of 
485.00°K  and  with  sudden  evaporation  of  2.5%  of  water  at  entry 
to 

(1)  stage  1 (curve  denoted  1) 

(2)  stage  2 (curve  denoted  2) 

(3)  stage  3 (curve  denoted  3) 

(4)  stage  4 (curve  denoted  4) 

(5)  stage  5 (curve  denoted  5) 

(6)  stage  6 (curve  denoted  6) 

In  examining  the  figures  it  is  clear  that  the  first  three  stages 
seem  to  fall  generally  into  one  category  and  the  latter  three  stages 
into  a second  category.  In  fact  the  characteristics  of  the  N.G.T.E. 
compressor  are  such  that,  as  shown  in  Fig.  2.3,  it  is  not  possible  to 
estimate  the  overall  performance  of  the  compressor  for  air-water  vapor 
operation  under  the  given  temperature  conditions  if  there  is  sudden 
evaporation  in  stages  1,  2 or  3. 

Several  observations  may  be  made  from  those  results. 

(1)  At  the  higher  rotational  speeds,  evaporation  of  2.5%  by  weight 
of  water,  when  it  occurs  at  entry  to  the  first  three  stages,  permits  only 
a restricted  range  of  operation  in  terms  of  corrected  mass  flow  at 
entry  to  compressor. 

(2)  At  the  higher  rotational  speeds,  evaporation  of  2.5%  by  weight 
of  water,  when  it  occurs  at  entry  to  stages  4,  5 and  6,  produces  a change 
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in  compressor  performance  but  there  is  still  a wide  range  of  mass  flows 
at  which  the  compressor  can  be  operated. 

(3)  The  differences  in  performance  that  are  obtained  when  water  is 
taken  to  undergo  evaporation  at  different  stations  in  the  compressor  are 
a clear  indication  of  the  influence  of  stage  loading. 

In  summary,  therefore,  the  effects  of  operating  a compressor  under 
the  condition  of  evaporation  of  water  at  some  station  in  the  compressor 
are  dependent  upon  stage  characteristics  (drop  in  efficiency  at  low  mass 
flows,  choking  at  high  mass  flows  and  the  mass  flow  range)  and  stage 
loadings . 

2. 1.3. 2.  The  results  of  calculations  on  the  Allison  Research  Compressor 
are  presented  in  figs.  2.4  - 2.8  for  operation  at  one  speed  of  rotation 
and  one  inlet  temperature  of  533. 2R. 

In  each  of  the  figures,  the  corrected  mass  flow  is  defined  as 
follows. 

corrected  mass  flow  = 

0 

■ */Vl7'W<Pol/Pref> 

where  the  temperature  has  been  normalized  with  respect  to  518. 4R  and  the 
pressure  has  been  normalized  with  respect  to  1.0  atmosphere  (standard). 
The  corrected  speed  is  defined  as  follows. 

N N 

corrected  speed  = — = 

* 

The  results  of  calculations  on  the  Allison  Research  Compressor 
are  presented  in  Figs.  2.4  - 2.8.  In  Fig.  2.4  the  individual  stage 
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characteristics  for  the  six  stages  of  the  compressor  are  presented  at 
one  speed  and  at  air  inlet  temperature  of  533. 2°R.  These  stage  char- 
acteristics have  been  employed  in  all  of  the  subsequent  calculations. 

In  Fig.  2.5,  the  overall  pressure  ratio  of  the  Allison  Research 
Compressor  is  presented  for  the  following  cases. 

(a)  Operation  of  the  compressor  with  air  only. 

(b)  Operation  of  the  compressor  with  sudden  evaporation  of  2.5 % 
of  water  at  the  stage  noted. 

In  Fig.  2.6  the  overall  temperature  ratio  is  given  for  the  same 
conditions  as  above. 

In  Fig.  2.7  the  overall  adiabatic  efficiency  is  given  for  the  same 
conditions  as  above. 

It  will  be  recalled  that  only  "natural  evaporation  cases"  have  been 
examined  here  and  also  at  one  speed  only.  Thus,  under  the  conditions  of 
operation,  there  is  no  possibility  of  sudden  evaporation  of  2.5%  of 
water  in  stages  1,  2 or  3.  The  first  location  where  such  evaporation 
can  occur  is  at  the  exit  of  the  third  stage  or  the  entry  to  the  fourth 
stage. 

The  influence  of  stage  loading  and  of  the  location  at  which 
evaporation  sets  in  is  clear  from  those  figures. 

In  Fig.  2.8,  in  order  to  illustrate  this  point  further,  the  per- 
formance characteristics  of  the  fourth,  fifth  and  sixth  stages  are 
presented,  each  stage  characteristics  corresponding  to  evaporation  at  the 
specified  stage.  Thus,  the  sixth  stage  performance  is  shown  when  sudden 
evaporation  occurs  at  entry  to  the  4th,  5th  and  6th  stages;  and  so  on. 
This  shows  the  cumulative  effect  of  various  downstream  stages  when 
evaporation  occurs  at  one  of  the  earlier  stages. 
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Thus,  in  Fig.  2.8,  the  following  notation  is  introduced. 

Curves  4/A,  5/A  and  6/A:  Performance  curves  of  4th,  5th  and  6th  stages 
with  100%  air. 

Curves  4/4,  5/5  and  6/6:  Performance  curves  of  4th,  5th  and  6th  stages 
with  100%  air  and  2.5%  steam  when  there  is  sudden  evaporation  at 
entrance  to  the  4th,  5th  and  6th  stages,  respectively. 

Curve  5/4:  Performance  of  5th  stage  with  100%  air  and  2.5%  steam  when 
there  is  sudden  evaporation  at  entrance  to  4th  stage.  This  curve, 
therefore,  represents  the  performance  of  the  5th  stage  when  both 
4th  and  5th  stages  operate  with  100%  air  and  2.5%  steam. 

Curve  6/4:  Performance  of  6th  stage  with  100%  air  and  2.5%  steam  when 
there  is  sudden  evaporation  at  entrance  to  4th  stage.  This  curve, 
therefore,  represents  the  performance  of  the  6th  stage  when  stages 
4,  5 and  6 operate  with  100%  air  and  2.5%  steam. 

Curve  6/5:  Performance  of  the  6th  stage  with  100%  air  and  2.5%  steam 
when  there  is  sudden  evaporation  at  entrance  to  5th  stage.  This 
curve,  therefore,  represents  the  performance  of  the  6th  stage  when 
stages  5 and  6 operate  with  100%  air  and  2.5%  steam. 

It  will  be  observed  that  both  5/4  and  6/4  are  closer  to  5/A  and 

6/A,  respectively,  than  5/5  and  6/6.  This  may  be  taken  as  an  indication 

of  the  cumulative  effect  of  stages  4 and  5 in  one  case  and  of  stages  4, 

5 and  6 in  the  other  case.  It  just  happens  in  this  particular  case  that 
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the  cumulative  effect  is  "favorable"  in  that  the  effect  of  addition  of 

4 

steam  has  decreased  as  the  fluid  has  gone  through  additional  downstream 
stages.  But  this  must  not  be  assumed  to  occur  universally.  In  another 
case  with  different  stage  characteristics  and  stage  loading,  an  entirely 
different  phenomenon  may  be  observed. 

2.2  Influence  of  Bleeding  or  Injection 

A simple  method  of  correcting  the  performance  of  a compressor  is 
through  arranging  bleed  or  injection  at  a chosen  stage  of  the  compressor.  \ 

Obviously  the  effect  of  bleed  or  injection  is  felt  in  the  stages  down- 
stream of  the  section  where  such  a correction  is  applied. 

We  now  consider  a multistage  compressor  operating  with  two  phase 
flow  and  suppose  the  evaporation  of  water  occurs  at  a certain  stage. 

If  the  performance  of  the  succeeding  stage  under  those  conditions  can  be 
brought  back  to  the  performance  that  the  compressor  would  have  with  air 
flow  only,  through  injection  or  bleeding,  then  the  succeeding  stages 
would  perform  as  with  air  alone  and  the  performance  of  the  compressor 
would  become  the  same  as  with  air  flow  under  those  conditions. 

We  consider  here  only  those  cases  where  the  water  content  by  weight 
and  hence  by  volume  is  quite  small.  Also,  we  restrict  attention  to 
steady  state  operation  sufficiently  far  from  stalling  conditions. 

When  bleeding  is  arranged  it  is  clear  that  the  blow-off  port  needs 
to  be  located  at  some  stage  downstream  of  the  stage  where  the  ingested 
water  is  likely  to  undergo  evaporation. 

When  injection  is  arranged,  one  has  to  decide  what  air  to  inject. 

One  of  the  options  --  and  this  may  be  the  only  practical  option  --  is 
to  divert  a small  amount  of  air  from  the  exit  plane  of  the  compressor 
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into  the  port  for  injection.  In  fact  this  is  the  case  considered  in 
the  following  analysis. 


2.2.1.  Problem  Formulation 

The  following  assumptions  are  introduced: 

(1)  Same  as  assumptions  (1)  under  2.1.1. 

(2)  The  air  that  is  injected  is  the  air  from  the  exit  plane  of 
the  compressor  when  it  is  operating  with  evaporated  water. 

Based  on  the  foregoing  assumptions,  the  problem  is  formulated  as 
follows.  Referring  to  Fig.  2.9,  we  consider  the  performance  of  a 
certain  stage  of  the  compressor  in  terms  of  corrected  mass  flow  at 
entry  to  the  stage  and  pressure  ratio  across  the  stage  at  a certain 
value  of  corrected  speed.  Suppose  the  performance  of  the  stage  with 
air  alone  is  designated  by  (A)  and  the  performance  of  the  stage  with 
evaporation  of  a certain  amount  of  water  is  designated  (A  + S).  It  is 
clear  that  there  is  a change  in  the  corrected  speed  of  the  stage 
following  evaporation  in  view  of  the  reduction  in  temperature  that 
arises  during  evaporation.  Suppose  now  that  the  particular  stage  is 
operating  at  point  (a)  on  (A)  with  air  alone  and,  on  evaporation  of  a 
certain  amount  of  water  at  entry  to  the  stage,  the  performance  of  the 
stage  shifts  to  point  (1).  The  problem  then  is  to  establish  the  bleed 
and/or  injection  (the  latter  with  the  mixture  drawn  from  the  exit  of 
the  compressor,  say)  that  is  required  in  order  to  readjust  the  perfor- 
mance of  the  stage  towards  point  (a  ).  However  point  (a  ) represents 
a certain  corrected  mass  flow  and  pressure  ratio.  In  practice,  in 
general,  it  may  be  impossible  to  return  the  compressor  stage  that  has 
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begun  to  operate  at  point  (1)  to  operate  at  point  (a).  In  a given  case, 
therefore,  it  is  appropriate  to  establish  whether  bleed  and/or  injec- 
tion has  to  be  introduced  at  (1  ) to  move  towards  (a)  and  what  the  effect 
of  the  magnitude  of  bleed  and/or  injection  is  at  (1). 

For  the  case  illustrated  in  Fig.  2.9  it  is  clear  that  bleeding 
alone  is  not  effective  in  changing  the  operating  point  (1)  towards  the 
operating  point  (a).  There  are  then  two  other  possibilities  for  adjust- 
ing the  operating  point  of  the  stage,  namely  injection  alone  or  a 
combination  of  bleeding  and  injection.  The  latter  may  become  necessary 
if  for  various  practical  reasons  one  sets  limits  on  the  magnitude  of 
injection  and  bleeding. 

2.2.2.  Solution  Procedure 

The  solution  procedure  and  the  associated  computer  program  (Purdue 
University  Compressor  Performance  Program,  PUCPP-002)  are  described  in 
Appendix  3 to  this  report. 

The  program  is  intended  for  recalculating  the  performance  of  a 
stage  with  bleeding,  injection  (or  air  at  any  desired  temperature  and 
pressure)  or  combinations  of  those. 

One  can  estimate  the  nature  of  the  change  in  mass  flow  required  for 
a desired  change  in  performance  or  the  change  in  performance  that  can  be 
realized  with  a given  change  in  mass  flow.  The  effect  of  injection/ 
bleeding  is  dependent  upon  the  stage  loading  in  addition  to  change  in 
mass  flow. 

2.2.3.  Results  of  an  Illustrative  Case 

Investigations  have  been  performed  on  the  4th  stage  of  the  N.G.T.E. 
#109  compressor. 
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Assuming  air  as  the  working  fluid,  the  performance  of  the  4th 
stage  of  the  compressor  is  established  for  a range  of  mass  flows  at 
entry  to  the  first  stage  of  the  compressor.  This  is  illustrated  in 
Fig.  2.10. 

Next,  assuming  that  2.5%  by  weight  of  water  has  evaporated  at 
entry  to  the  4th  stage  the  performance  of  the  4th  stage  is  established 
again  and  is  also  shown  in  Fig.  2.10. 

It  may  be  pointed  out  here  that  the  performance  curves  shown  in 
Fig.  2.10  for  stage  4 are  a consequence  of  the  range  of  mass  flow 
selected  at  entry  to  the  1st  stage  of  the  compressor.  The  available 
stage  characteristics  for  the  1st,  2nd  and  3rd  stages,  even  after  some 
extrapolation,  do  not  permit  any  larger  mass  flows  to  be  considered, 
than  used  here,  at  entry  to  the  1st  stage. 

Now,  calculations  were  performed  to  establish  the  effect  of 
injection  of  three  different  mass  flows  from  the  exit  of  the  compressor 
to  the  entry  section  of  the  4th  stage  when  that  stage  was  operating 
with  2.5%  water  evaporated  at  entry  to  it.  The  magnitudes  of  the 
injected  mass  were  2.0%,  5%  and  10%.  The  resulting  performance  of  the 
4th  stage  is  again  illustrated  in  Fig.  2.10. 

Two  considerations  to  be  taken  into  account  in  performing  such 
calculations  are  as  follows. 

(1)  The  effect  of  injecting  a certain  percent  of  compressor  exit 
flow  into  an  upstream  stage  such  as  the  4th  stage  can  be  established 
only  by  iteration  since  the  compressor  exit  condition  itself  changes 
with  upstream  injection. 
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(2)  From  the  same  reasoning,  the  compressor  exit  condition 
changes  when  different  amounts  of  exit  flow  are  introduced  into  an 


upstream  stage  such  as  the  4th  stage. 

Referring  to  the  operating  lines  presented  in  Fig.  2.10,  while  it 
appears  that,  with  5%  injection  of  compressor  exit  flow,  it  is  possible 
to  come  close  to  the  performance  characteristic  with  air  alone,  there 
are  large  differences  between  the  individual  operating  points  with 
air  alone,  with  evaporation  but  no  injection  and  with  evaporation  and 
injection. 

It  is  important  to  observe  here  that  the  foregoing  applies  to  a 
particular  stage  in  a given  compressor  with  a chosen  magnitude  of 
evaporation  at  entry  to  the  stage.  The  same  conclusions  may  not  be 
reached  in  another  case.  All  that  may  be  said  at  this  stage  is  that, 
there  are  limits  to  the  magnitude  of  bleeding  and/or  injection,  the 
stage  loading  and  the  location  and  magnitude  of  evaporation  will 
determine  the  nature  and  extent  of  the  changes  that  can  arise  through 
bleeding  and/or  injection. 

2.3  Influence  of  Location  of  Hater  Vaporization 

In  both  of  the  previous  discussions  under  Sections  2.1  and  2.2, 
we  have  emphasized  the  location  of  water  vaporization  in  addition  to 
the  mass  of  water  evaporated  and  the  stage  loading  of  the  stages 
following  the  location  of  vaporization. 

In  Ref.  3,  a test  has  been  reported  wherein  a slug  of  water  was 
ingested  into  a compressor  under  two  different  conditions:  (l)  cold 
condition  and  (2  ) post  warm-up  condition.  It  was  found  that  under 
cold  conditions  there  was  no  noticeable  effect  on  the  engine  performance. 
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However,  when  the  slug  of  water  was  ingested  after  warm-up  of  the  engine, 

4 

there  was  both  surging  of  the  compressor  and  after-burner  blow-out. 

The  warm-up  condition  may  be  considered  as  equivalent  to  a rise  in 
the  temperature  of  ambient  air  or  as  equivalent  to  a shift  upstream  of 
the  location  of  water  vaporization  in  a given  compressor. 

In  the  alternative,  the  warm-up  condition  may  have  lead  to  instan- 
taneous differential  expansion  between  the  rotor  and  the  casing  of  the 
compressor  and  hence  to  enlarged  clearance  and  consequent  change  in 
performance.  At  this  stage,  it  is  not  possible  to  state  positively 
what  the  cause  for  the  changed  performance  is.  Nevertheless,  the  per- 
formance of  a compressor  has  been  shown  here  to  deteriorate  when 
evaporation  of  water  is  advanced  upstream  in  the  compressor. 
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Figure  2 2 Pressure  Ratio  vs.  Mass  Flow — 8000  RPM 


ALLISON  RESEARCH  COMPRESSOR 


Figure  2.5  Pressure  Ratio  vs.  Corrected  Mass  Flow— 9939  RPM 

(N//S‘=9800) 


ALLISON  RESEARCH  COMPRESSOR 


Temperature  Ratio  vs.  Moss  Flow— 9939  RPM  (N/V@  = 9800) 


ALLISON  RESEARCH  COMPRESSOR 


Adiabatic  Efficiency  vs.  Corrected  Mass  Flow  — 9939  RPM  (NA/S=9800) 


ALLISON  RESEARCH  COMPRESSOR 
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pdrating  Condition 


SECTION  III 


FURTHER  DEVELOPMENT  OF  TWO-PHASE  MODEL 

It  will  be  recalled  that  in  Ref.  1 the  presence  of  water  was  ac- 
counted for  in  terms  of  a continuum  of  droplets,  but  with  the  assump- 
tion that  the  non-random  distribution  of  droplets  does  not  necessitate 
the  inclusion  of  a separate  contribution  of  pressure  due  to  the  droplets. 
The  droplets  were  assigned  an  averaged  size  with  averaged  properties. 

It  is  now  proposed  to  improve  that  model  in  two  respects. 

(1)  Identification  of  three  droplet  sizes;  and 

(2)  Accounting  for  interchange  of  population  densities  of 
droplets  between  the  three  classes. 

In  this  fashion  it  is  hoped  to  account  for  the  presence  of  a variety 
of  discrete  droplets  by  means  of  three  classes  of  droplets,  each  asso- 
ciated with  its  own  averaged  size  and  averaged  characteristic  properties. 
The  choice  of  three  classes  (rather  than  some  larger  number)  is  arbitrary 
but  may  prove  adequate,  especially  in  view  of  other  inherent  uncertainties. 

3.1 . Droplet  Classes 

A small  study  has  been  performed  to  determine  the  droplet  sizes 
(range  and  distribution)  that  might  be  considered  as  most  representative 
during  rain  ingestion.  This  study  is  based  on  Refs.  4-6  • 

The  size  ranges  of  atmospheric  water  droplets  is  shown  in  Fig.  3.1  . 

The  size  distributions  of  water  particles  in  rain  and  clouds  is  shown 
in  Figs.  3.2  and  3.3,  respectively.  The  values  indicated  in  those 
figures  should  be  considered  as  typical,  noting  that  there  can  be  wide 
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discrepancies  in  the  estimates  of  such  distributions.  Utilizing  the 
information  in  Figs.  3.2  and  3.3,  a particular  classification  of 
water  particles  has  been  chosen,  as  shown  in  Fig.  3.4,  for  use  in  this 
report. 

It  will  be  observed  that  three  ranges  of  particle  sizes  are  selected 
as  follows. 

(1)  drizzle,  ranging  from  5 to  30  microns; 

(2)  small  droplets,  designated  Class  1 and  ranging  from  100 
to  800  microns;  and 

(3)  large  droplets,  designated  Class  2. 

Those  three  classes  are  further  characterized  as  follows. 

(a)  In  the  case  of  drizzle,  it  is  assumed  that  it  is  a continuum 
moving  with  the  air  with  no  viscous  interaction  identified  between  the 
particles  and  the  air.  These  small  particles  are  further  assumed  to 
introduce  no  pressure  on  the  system. 

(b)  In  the  case  of  the  other  two  sizes  of  particles  it  is  assumed 
that  each  class  may  again  be  treated  as  a continuum  of  some  assigned, 
averaged  size  particles  but  each  class  moving  with  its  own  characteristic 
velocity.  In  view  of  the  larger  sizes  of  these  particles,  viscous  in- 
teraction is  taken  into  account  between  the  particles  and  air.  Again 

it  is  assumed  that  neither  of  these  classes  of  particles  exerts  any 
pressure  force  on  the  system. 

(c)  Finally,  in  the  case  of  the  two  larger  size  particles,  it  is 
assumed  that  coalescence  and  disintegration  of  droplets  can  take  place 
according  to  accepted  criteria  for  such  processes.  However,  the  small 
droplets.  Class  1,  when  they  coalesce  (among  themselves  or  with  Class  2 
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droplets)  are  always  expected  to  yield  Class  2 droplets.  Similarly  the 
large  droplets.  Class  2,  when  they  break-up  are  always  expected  to  yield 
Class  1 droplets.  In  other  words.  Class  1 droplets  cannot  coalesce  into 
other  Class  1 droplets,  nor  can  Class  2 droplets  break-up  into  other 
Class  2 droplets. 

3.2.  Model  for  Coalescence  and  Break  Up 

The  interaction  between  Class  1 and  Class  2 droplets  is  expected 
to  arise  either  through  droplet  coalescence  or  droplet  break  up.  The 
interaction  model  is  similar  to  that  employed  in  elementary  molecular 
theory  (Ref.  7). 

The  model  is  developed  as  follows. 

(a)  The  droplets  in  Classes  1 and  2 each  lie  between  the  diameters 
Di  min  and  Di  max  ^ = ^»2)»  but  c^ass  can  be  represented  by  a character- 
istic size  D. . 

(b)  Suppose  N.j  is  the  number  of  droplets  of  range  i per  unit  volume. 
A distribution  function  f..  is  then  introduced  such  that 

dNi  = N.f.dD  (1) 

where  dNi  is  the  number  of  droplets  of  type  i per  unit  volume  in  the 
diameter  range  and  D..  + dD. 

(c)  The  velocity  of  all  of  the  droplets  of  type  i is  denoted  by 

V 

(d)  Coalescence  is  supposed  to  occur  by  collisions  between  droplets 
of  the  two  Classes  1 and  2,  but  not  by  collisions  in  Class  1 or  2. 

Break  up  is  assumed  to  occur  when  a critical  value  of  Weber  number  is 
reached . 
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3.2.1.  Coalescence  model 

A notion  of  a collision  time  At  is  introduced,  for  example  At  can 
be  related  to  the  half  period  of  a free  oscillation  of  a spherical  drop- 
let. Then  the  collision  volume  or  the  volume  traversed  by  droplets  during 
a collision  can  be  expressed  by  the  relation, 

AV  = J (D}  + D2)2(Ul  - u2)At.  (?) 

Let  n.j  be  the  number  of  collisions  experienced  by  particles  of 
Class  2,  which  is  the  same  as  that  experienced  by  particles  of  Class  1, 
in  the  volume  AV.  Then,  it  follows  that 

^ = N] At  | } (D-,  + D2)2(u1  - u2)f]dD  ( 3) 

We  will  assume  here  that  the  mean  droplet  size,  the  mean  and  standard 
deviations  are  such  that,  when  the  droplets  obey  normal  distribution, 
nearly  all  of  the  collisions  will  have  been  accounted  for  in  the  analysis. 

Now,  the  total  number  of  collisions  per  unit  time  can  be  written  as 

n = Nj/At  for  ^2  > a7  (4) 

N N2  1 

= At  * (1/AV)  f0r  ,N2  - AV  ( 

or,  in  general,  as 


n = RiN2/At  where  R = N-j/N2  for  N2  > 1/AV  ( £.) 

R = N-|  - AV  for  N2  <.  1/AV 

We  can  then  express  the  loss  rate  of  mass  of  type  1 droplets  per 
unit  volume  of  mixture  due  to  collisions  in  the  following  form. 


RN2p-|TT 


D3f i dD 


(7) 
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It  may  be  pointed  out  that  the  loss  rate  of  type  1 droplets  is  equal 
to  the  gain  rate  of  type  2 droplets. 

Let  ai  be  the  bulk  volume  concentration  of  type  i droplets  in  the 
mixture.  Then 

ai  = l ^i  \ Di  fi  dD  (8> 

The  mass  density  of  droplets  of  type  i can  then  be  written  as  follows. 


p 


i 


aip£ 


(9) 


where  is  the  density  of  the  liquid. 

Finally,  it  is  necessary  to  identify  the  number  of  droplets  of  each 
class,  say  m^ , that  would  become  created  (due  to  condensation)  or  loss 
(due  to  evaporation)  per  unit  volume  of  mixture.  We  can  define  the  mass 
of  such  droplets  as  follows. 


"Gi  = (l"ci  ' "ei1 


(10) 


where  the  subscripts  c and  e refer  to  condensation  and  evaporation  re- 
spectively. 

On  a continuum  basis,  therefore,  one  has  to  examine  the  balance 

between  p.,  m„.  and  m_..  This  however  does  not  take  into  account  the 
1 LI  ul 

break  up  of  droplets,  which  will  be  discussed  in  the  following  section. 


3.2.2.  Break-up  model 

In  order  to  account  for  droplet  break-up,  we  introduce  the  notion 
of  Weber  number,  We,  defined  by 

We  = (11) 
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where  p and  a are  the  environmantal  pressure  acting  on  the  droplet  and 
surface  tension,  respectively. 


We  assume  that  only  droplets  of  Class  2 can  undergo  break  up. 

The  environmental  pressure  acting  on  the  droplet  can  be  written  as 
the  sum  of  aerodynamic  force  and  collisional  force, 

P 


aero  ' ■>,  I 


and 


where 


, coll 

3coll  - A* 


F = tv  (u,  - u0)n 
col  1 At  1 1 2'  'c 


(12) 


(13) 


(14) 


and  m-j  is  the  average  mass  of  a droplet  of  Class  1 and  nc  is  the 

★ 

"efficiency"  of  collisions.  The  area  normal  to  F ^ is  given  by  A . 

Now,  suppose  that  the  number  of  droplets  breaking  up  per  unit  volume 
is  Nb-  Then,  introducing  an  arbitrary  constant  B, 


N 


b = A exp  (B  We) 


(15) 


One  can  adopt  the  following  boundary  conditions  for  Eq.  4.12,  namely 


^b  = ^*bo  at  = We  (m''n''mnum) 


and 


Nb  = 1 at  We  = We  (critical) 
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Thus,  it  follows  that 


Nb0  + exp  [We/(Wemin  - Wecr)] 
b - Nb0  * exp  LWecr/(Wemin  - Hecr)J 


(16) 


The  total  number  of  droplets  generated  per  unit  volume  per  unit 
time  due  to  break  up  is  given  by 


N 


lb 


= JL 

AT 


N.VdN 

b 


(17) 


where  K = 0,1  for  We  < We^  and  We  _>  Wecr-  Thus,  the  gain  of  mass  of 
type  1 droplets  due  to  a break  up  of  type  2 droplets  is  given  by  N^rn^. 
This  of  course  is  equal  to  the  loss  of  mass  of  type  2 droplets  due  to 
break  up. 

We  can  now  conclude  that  there  should  be  a balance  between  , ri^.»  mCl- 

and  mBi  taking  account  of  condensation  and  evaporation,  coalescence  and 
break  up. 


3.3.  Mass,  Momentum  and  Energyvof  Liquid  Phase 

As  stated  earlier,  the  part  of  the  liquid  phase  that  is  in  the  form 
of  drizzle  is  assumed  to  behave  as  part  of  the  gas  phase.  It  will  under- 
go evaporation  when  conditions  are  right  but  there  is  no  viscous  inter- 
action between  the  drizzle  particles  and  air,  and  no  interaction  among 
the  drizzle  droplets. 

In  essence,  therefore,  only  the  Class  1 and  Class  2 liquid  particles 
enter  into  the  discussion  of  the  liquid  phase. 
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3.3.1,  Mass  balance  equation 


The  mass  balance  equation  for  Class  i droplets  may  be  set  up  as 
follows. 


8pi  a 


m, 


'Gi 


+ mCi  + m' 


Bi 


(13) 


3.3.2.  Momentum  balance  equation 

In  addition  to  the  forces  taken  into  account  in  setting  up  the  momen- 
tum balance  equation  for  the  liquid  phase  in  Ref.  1,  we  have  to  take  into 
account  the  following  forces  when  two  classes  of  droplets  are  introduced. 

(1)  momentum  change  due  to  mass  change  of  each  class  of  droplets; 
and 

(2)  momentum  change  because  of  interaction  between  the  two  classes 
of  droplets. 

The  latter  arizes  because  of  each  class  of  droplets  being  treated 
as  a continuum  and  the  interaction  between  the  two  droplet  continua. 

We  shall  denote  this  by 

3.2.3.  Energy  balance  equation 

In  addition  to  the  energy  components  taken  into  account  in  Ref.  1, 
we  have  to  include  the  following  energy  components. 

(1)  Energy  change  due  to  mass  change  of  each  class  of  droplets; 

(2)  Energy  required  in  connection  with  forces  and 

(3)  Energy  required  for  the  change  of  droplet  surface  area  during 
formation  of  new  droplets. 


Regarding  (iii),  if  E^ ; is  the  energy  involved,  we  may  write  this 
as 

Esi  = £ [ 7rD2f2dD]  (19) 

3.4#  Equations  of  Motion 

The  equations  of  motion  are  presented  in  Appendix  . The  manner  in 
which  they  are  presented  is  the  same  as  in  Ref.  1,  namely, 

I.  Three  dimensional  equations  of  motion  in  intrinsic  coordinates. 

II.  Axisymmetric  equations  of  motion  in  intrinsic  coordinates. 

III.  Three  dimensional  equations  of  motion  in  cylindrical  coordinates. 

IV.  Axisymmetric  equations  of  motion  in  cylindrical  coordinates. 
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SECTION  IV 


SOME  MODIFICATIONS  TO  THE 
COMPRESSOR  PROGRAM  UDO-300 

The  original  program  (Ref.  8)  has  been  developed  for  air  flow  through 
a compressor.  In  order  to  make  it  applicable  to  air-water  mixture  flow, 
it  is  essential  to  incorporate  various  modifications.  Two  of  such  modi- 
fications have  been  developed  to-date  and  they  pertain  to  the  following. 

(1)  Need  for  taking  into  account  variable  specific  heats  and  mole- 
cular weight  along  the  flow  path  for  the  gas  phase  of  the  mix- 
ture; and 

(2)  Need  for  taking  into  account  variable  acoustic  velocity  for  j 

the  air-water  mixture  along  the  flow  path. 

4.1  Variable  Specific  H eats  and  Molecular  Weight 
The  gas  phase  in  the  compressor  is  assumed  to  be  a mixture  of  per- 
fect gases.  Any  water  vapor  generated  due  to  evaporation  of  water  is 
assumed  to  mix  instantaneously  and  homogeneously  with  the  surrounding 
air. 

The  specific  heat  at  constant  pressure  for  the  gas  phase  is  thus  a 
function  of  the  local  temperature. 

The  molecular  weight  for  the  gas  phase  is  a function  of  the  mass 
fractions  of  air  and  water  vapor  (or  steam). 

The  changes  required  in  the  FUNCTION  subprograms  of  UDO-300  are 
described  in  Appendix  5. 

4 . 2 Acoustic  Velocity  in  2-phase  Flow 
As  stated  in  Ref.  1,  acoustic  velocity  in  a two-phase  mixture  (with 

one  discrete  phase)  can  be  defined  as  follows. 
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/I  - a \ a i/2 

Cm  = l{{]  ' °v)pg  + °vpp}  * {(  VV  + (20) 

pgcg  ppcp 

where  cm  is  the  mixture  sonic  velocity,  a , the  (particulate)  liquid 
volume  fraction,  p,  the  density  and  ( ) and  ( )p  refer  to  the  gas  and 
liquid  phase  respectively.  That  relation  does  not  include  the  effect 
of  size  distribution  of  water  drops.  In  order  to  account  for  that  and 
also  the  nonlinear  interactions  between  attenuation  and  scatter  of  sound, 
empirical  corrections  can  be  introduced. 

Details  of  the  modifications  to  the  FUNCTION  subprograms  are  given 
in  Appendix  5. 
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SECTION  V 
DISCUSSION 


The  Report  describes  five  aspects  of  investigations  in  the  subject 
of  water  ingestion  in  axial  flow  compressors. 

(1)  Analysis  of  performance  of  selected  compressors  with  water 
ingestion. 

(2)  Determination  of  the  influence  of  blow-off  and  injection  into 
a compressor  operating  with  air-water  mixture. 

(3)  Development  of  a model  for  an  engine  operating  with  water 
ingestion  into  the  compressor. 

(4)  Formulation  of  flow  equations  for  an  axial  flow  compressor 
operating  with  air-water  mixture  accounting  for  several  classes 
of  droplets  and  such  processes  as  droplet  coalescence,  break- 
up and  evaporation. 

(5)  Some  simple  modifications  to  the  UD0-300  program  to  account 
for  variable  molecular  weight,  specific  heats  and  acoustic 
velocity. 

5.1  Conclusions 

(1)  The  precise  effects  of  the  presence  of  water  and  of  its  evapor- 
ation on  compressor  characteristics  and  surging  depend  upon 
(a)  the  stage  loading,  (b)  the  change  in  the  performance  due 
to  presence  of  water,  (c)  the  location  of  evaporation,  and 
(d)  the  changes  introduced  in  the  compressor,  locally  and 
upstream,  due  to  the  various  processes  associated  with  water 
ingestion  and  evaporation.  In  an  approximate  analysis  where 
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(a)  the  presence  of  water  in  liquid  form  has  been  assumed  to 
have  little  effect  on  compressor  performance,  (b)  a finite 
quantity  of  water  is  expected  to  undergo  sudden  evaporation 
at  predetermined  stages  and  (c)  upstream  effects  induced  by 
changes  further  downstream  are  neglected,  it  has  been  shown 
that  (a)  the  location  of  evaporation  has  a determining  effect 
on  the  deterioration  of  the  overall  performance  of  the 
compressor  and  (b)  such  deterioration  of  the  compressor  depends 
upon  the  stage  loading  and  characteristics  of  the  compressor. 

(2)  In  general,  it  is  extremely  difficult  in  practice  to  correct 
the  performance  of  a compressor  that  is  operating  with  water 
evaporated  in  it  (at  some  intermediate  stage)  by  means  of 
blow-off  or  injection.  In  selected  compressors,  it  has  been 
found  that  the  effects  of  evaporation  of  even  2.5%  by  weight 
of  water  cannot  be  rectified  by  injection  of  5%  by  weight  of 
air  flow  at  practicable  temperature  and  pressure. 

5.2  Possible  Directions  for  Further  Studies 

(1)  There  is  scope  for  improving  the  one-dimensional  analysis 
currently  employed  to  assess  the  performance  of  selected 
compressors  operating  with  air-water  mixture  flow. 

(2)  The  modeling  of  overall  engine  performance  under  nonsteady 
state  conditions  when  there  is  water  ingestion  is  of  crucial 
importance  in  understanding  the  effects  of  water  ingestion  on 
control s. 

(3)  The  flow  equations  need  to  be  reformulated  in  the  regions  of 
blade  and  casing  surfaces  and  also  to  account  for  mixture 
compressibility  effects. 
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(4)  Most  important  of  such  further  studies  is  the  one  that  will 
provide  proper  and  adequate  experimental  information  on  the 


effects  of  water  ingestion  into  isolated  compressors  and 
engines.-  Instrumentation  and  the  development  of  observation 
techniques  for  such  studies  require  in  most  cases  further 
development . 


APPENDIX  A 

PERFORMANCE  PARAMETERS 

The  performance  of  a compressor  may  be  expressed  in  terms  of  the 
following  parameters  on  a one-dimensional  flow  basis. 

1)  Va/a 

2)  U/a  and 

3>  VP0i 

4)  nc 

where  V,  a,  1),  and  PQ  and  nc  represent  the  axial  velocity,  acoustic 
speed,  rotational  speed,  total  pressure  and  efficiency,  and  ( )i  and 
( ) refer  to  the  inlet  and  exit  conditions  of  a compressor.  If  a 

single  stage  of  a multi-stage  compressor  is  considered,  one  can  again 
express  the  performance  on  the  basis  of  the  same  parameters  for  the 
stage  under  consideration. 

In  general,  the  parameter  V /a  may  be  expressed  in  terms  of  the 

a 

following  relation. 


Y + 1 


If  a compressor,  designed  originally  for  operation  with  air  (a),  is 
expected  to  be  operated  with  another  gas  (g),  it  follows  that 
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where 


u 


J- 


+ 1 


Y-I 


M2] 


2(y  - T) 


at  corresponding  operating  points. 

Similarly,  the  parameters  U/a  may  be  expressed  in  terms  of  the 
following  relation. 


M2] 


1/2 


(A-3 


Again,  considering  the  operation  of  the  same  compressor  with 
different  gases,  it  follows  that 


(A  -4 


where  [ ] = • [1  + Y-  y -1-  M2]1/2 


at  corresponding  operating  points. 

It  may  be  pointed  out  that  under  the  one-dimensional  flow  approxima- 
tion corresponding  operating  points  are  assumed  to  yield  the  same  pressure 
ratio  and  efficiency.  In  order  to  take  into  account  blade  passage, 
blade  clearance  and  inter  blade-row  spacing,  a more  sophisticated 
approach  than  one-dimensional  flow  is  required. 

It  is  usual  practice  to  express  the  and  parameters  as 


y-i: 
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corrected  quantities  by  using  reference  standard  values  of  inlet- 
pressure  and  temperature.  The  corrected  mass  flow  parameter  and 
rotational  speed  parameter  for  a compressor  are  designated  as 


and 


ir- 


respectively. 

Throughout  the  analysis,  unless  it  is  specified  that  the  perfor- 
mance of  a particular  stage  is  being  presented  with  respect  to  its  own 
stage  inlet  conditions,  all  of  the  compressor  performance  maps  are 
presented  in  terms  of  corrected  mass  flow  and  corrected  speed. 
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APPENDIX  B 


PURDUE  UNIVERSITY  COMPRESSOR  PERFORMANCE  PROGRAM-PUCPP-OOl 

The  purpose  of  this  program  is  to  calculate  the  stage  and  overall 
compressor  performance  using  a one-dimensional  technique.  The  analysis 
is  limited  to  a six-stage  compressor  for  which  exists  stage  data  relating 
axial  velocity  to  stage  pressure  rise  and  stage  efficiency.  The  following 
four  sections  describe  the  necessary  input  data,  the  resulting  output 
data,  the  one-dimensional  method  and  capabilities  of  the  program,  and 
the  program  listing. 

1.  Input  Data:  There  are  two  items  of  input:  (i)  stage  performance  and 
geometry  data,  and  (ii)  inlet  and  operating  conditions.  Both  integer  and 
real  input  formats  are  used.  Integer  numbers  are  punched  in  fields  of 
five  locations  starting  with  column  1.  No  decimal  points  may  be  used. 

Integer  fields  are  right-hand  justified.  Real  numbers  are  punched  in 
fields  of  ten  locations  starting  with  column  1.  Decimal  points  must  be 
used  for  real  numbers  if  right-hand  justification  in  each  field  is  to  be 
ignored. 

A listing  of  the  input  cards  is  given  in  the  following  chart,  where 
one  line  represents  one  card. 

N NDATA 

RR  A j 

Occurs  N times 

PHI  SI  ETA  - Occurs  NDATA  timesj 

NFLOW 

FLOW  RPM  T01  P01  CP  W1  - Occurs  NFLOW 

times 


Variable  Definitions  - The  input  data  must  have  units  corresponding 
to  the  dimensions  given  below  (where  applicable). 

Length  feet 

Time  seconds 

Mass  pound-mass 

Force  pound-force 

Temperature  Degrees  Rankine 

The  one  exception  to  this  is  the  input  value  of  the  compressor  rotational 
speed,  revolutions  per  minute. 

N Number  of  compressor  stages  N _<  6. 

NDATA  Number  of  data  points  to  be  input  that  give  stage 

performance  information.  There  must  be  NDATA  data 
points  for  each  stage  NDATA  _<  15. 

RR  Radius  at  which  calculation  is  to  be  made.  This  may 

be  an  average  radius  for  the  stage.  There  must  be  N 
values  of  RR  input.  Units: feet. 

A Stage  inlet  cross-sectional  area.  There  must  be  N 

values  of  A input.  Units:feet2. 

PHI  Stage  performance  information  - stage  inlet  axial 

velocity/rotor  speed  (determined  at  RR).  There  must 
be  NDATA  values  of  PHI  input  for  each  stage. 

SI  Stage  performance  information  - stage  pressure  rise/ 

stage  inlet  dynamic  head.  SI  must  correspond  to 
each  value  of  PHI . 

ETA  Stage  performance  information  - stage  adiabatic 

efficiency.  ETA  must  correspond  to  each  value  of 
PHI. 

It  is  necessary  to  input  the  stage  performance  information  (PHI, SI  and  ETA) 
in  ascending  order  beginning  with  the  smallest  value  of  PHI,  for  each  stage. 

NFLOW  Number  of  mass  flows  to  be  input. 

FLOW  Operating  mass  flow.  There  must  be  NFLOW  values  of 

FLOW  input.  Each  value  of  FLOW  must  have  the  remaining 
variables  associated  with  it.  Units:pound  mass/seconds. 
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RPM  Compressor  rotational  speed.  Must  correspond  to 

each  value  of  FLOW.  Units:minute_1 

T01  Compressor  inlet  total  temperature.  Must  correspond 

to  each  value  of  FLOW.  Units :Degrees  Rankine. 

P01  Compressor  inlet  total  pressure.  Must  correspond  to 

each  value  of  FLOW.  Units: pound-force/feet2. 

CP  Specific  heat  at  constant  pressure.  Assumed  to  be 

comstant  throughout  compressor.  Must  correspond  to 
each  value  of  FLOW.  Units:feet-pound-force/pound- 
mass/Deg.  Rankine 

W1  Compressor  inlet  molecular  weight.  Must  correspond 

to  each  value  of  FLOW.  Units :pound-mass/pound-mole 

2.  Output  Data:  There  are  two  types  of  output  data:  (i)  compressor 
performance,  and  (ii)  input  errors.  Providing  no  input  errors  have  occurred, 
the  program  will  give  values  for  the  following. 
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1.  Compressor  Total  Pressure  Ratio 

2.  Compressor  Total  Temperature  Ratio 

3.  Compressor  Adiabatic  Efficiency 

4.  Stage  Performance  Characteristics 

- stage  total  pressure  ratio 

- stage  total  temperature  ratio 

- stage  adiabatic  efficiency 


In  addition,  the  input  inlet  operating  conditions  corresponding  to  these 
performance  values  will  be  printed. 

In  the  event  that  the  input  inlet  operating  conditions  cause  the  stage 
performance  to  exceed  the  limits  of  information  (PHI,  SI  and  ETA)  that  has 
been  input  for  that  stage,  no  extrapolation  is  carried  out.  Instead,  the 
analysis  for  that  operating  point  ceases  and  the  program  continues  with 
the  next  operating  point. 
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The  following  message  is  printed: 


STAGE  X DATA  IS  NOT  SUFFICIENT  - CALCULATION 
TERMINATED  FOR  Y MASS  FLOW, 

where  X is  the  stage  where  the  limit  was  exceeded  and  Y is  the  mass  flow 
under  consideration.  If  the  compressor  performance  for  that  mass  flow 
is  desired,  the  stage  performance  information  for  stage  X must  be  expanded. 
3.  General  One-Dimensional  Method:  The  analysis  is  intiated  with  the 
input  values  of  the  inlet  operating  conditions.  A value  of  4>  is  deter- 
mined which  is  then  compared  with  the  stage  performance  data  input  to  the 
program  (PHI). 

^calc  va^u  pA2TrrRPm 


upon  determining  two  values  of  PHI  stored  in  input  files  that  flank  the 
calculated  value  of  't'caic»  a weighting  factor  is  determined: 


Wt.  Fact.  = (<t> 2 - <l>calc^^2  ' ‘h  ^ 


where  = immediately  greater  than  4>caic 
0 = PHI  immediately  less  than  <p  . 

Cd  I C 

It  is  assumed  that  a straight  line  connects  the  points  input  as  stage  per- 
formance information.  Thus,  with  the  above  weighting  factor  it  is  possible 
to  determine  the  values  of  ip  and  n (input  as  SI  and  ETA,  respectively) 
corresponding  to  <t>calc  Thus, 

^calc  = ( ^ ‘ VWt-  Fact-  + ^2  ( B-3 
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and 


ncalc  = ^nl  " n2^  Wt‘  Fact'  + n2’  (P-4 

where  (4^  and  t^)  and  (^2  and  n2)  correspond  to  (4>  1 ) and  (d> 2) , respectively. 
The  stage  pressure  rise  is  found  from  ^ca-|c: 

iPo  ' W'  (,/2pVa>  (B'5 


where  V is  the  stage  inlet  axial  velocity.  The  stage  temperature  ratio 

d 

is  found  by,  (i)  assuming  isentropic  flow, 


/T 


■ <y  v 


(y-1)/y 


(B-6 


(ii)  then  accounting  for  losses  by  the  use  of  the  stage  adiabatic  efficiency 
in  the  following  manner: 


where,  AT 


0, 


^actual  = 4T°isentropic/nad1ab 
is  found  using  equation  (6),  and  (iii) 


( B-7 


isentropic 


(t„  nr 


) = (T  + AT  )/T  . 


°2  °1  actual  U1  ^actual  U1 


( B-8 


The  stage  and  compressor  adiabatic  efficiency  is  determined  by, 

(P0  /P0  )(y_1)/y  - 1 
- °2  °1 

nadiab  T /T  - “1  ’ 

°2  °1 


where  2 represents  the  exit,  and  1 the  inlet,  of  the  stage  or  compressor, 
respectively. 
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1.  dimension  phko.is*.  sks.ism  etxia.ism  tout*.  aiat 

».  UIhLNSION  SPHlfcl.  STRIP).  SfTAlP) 

3.  R»1S4S,0 

«,  0 = 32.17** 

3.  Plt*3.1*lS»2T 

4.  WRITE  (4.1031 

T.  HERO  13.110)  N.nOATA 

6.  00  102  1 = 1. N 

9.  READ  13.1111  RRlD.Ad) 

10.  00  lol  HSl. NOMA 

11.  READ  (3.112)  PH]d.H).SKl*K>.tTAI|,M 

12.  101  CONTINUE 

13.  102  CONIINUC 

14.  READ  (3.113)  NFLOW 

15.  00  106  L=l. NFLOW 

14.  ST0H1M.0 

17,  ST0H2=1 .0 

16.  head  (3.11<*>  RPn.TOl.POl.cP.Wl .FLOW 

19.  OANAsl .0/1 1.0>N/ICP=W1 > ) 

20.  OanAF<(r.AH«-1.0)/r,AnA 

21.  DO  103  1=1, N 

22.  U:2.(i>('IC*HM(I|,RPH/40.0 

23.  Hh0=P01*W1/(H.T01I 

21.  VA=FLOW/IRHO*AI I > I 

23.  PHIKVA/U 

24,  K*1 

27,  IF  ( (PMlI-PHld,  1||, CE. 0.01  00  TO  104 

2ft.  103  WH1TL  (4.113)  I, FLOW 

29.  00  TO  106 

30.  104  KiH»l 

31.  ! F 'K.EO.INOATAd))  00  TO  103 

.'7.  1*  (PHI. 'hIU.ki.GT. 0.01  00  TO  104 

ST.  F ."'.Txf  Ml  I I ,*■)  -PHM  )/(CHT  ( I ,K)-PHI  (1  ,K-1)  | 

3".  Sll  = -F  TCT.  (0:  . 1 1 1 .K-l)  MSI  d.H) 

33.  ET»I  = -F*C  • * 1 L T ,1  I .Hi-LTAf  1 .K-l)  META  I I.  Hi 

34.  UELP  = S1 1 •0.3>Hmo«i.>  *2/0 

3,.  p02=pPi»0ELP 

IF,  4P((I i=302/POl 

39.  T R*  ( S*  H 1 1 ) ) ••GAHAr 

40.  0ELTa=(  T«»Tol-Tol I/lTAI 

41.  T02A=T01*OCLTA 

42.  Sf(i) * = 1 U2A/T01 

<3.  G£  i«(I  l = (=PHd  )..OAHAF-1.0l/(STR(lMl.0l 

44.  P01«P02 

43.  103  T01*T02A 

44.  00  104  1=1, N 

47.  TPHaSPp ( 1 ) aSIORl 

46.  TTH:STN(1)=ST0R2 

40.  STOH 1 =TPH 

30.  104  ST0«2  = T TR 

31.  UTA  = (TI‘H»=GAHAF-l,0)/tTTH-l.O) 

32.  WHITE  (4,11b)  FL0w.NPh.T01.P01.Wl 

33.  WH11L  (4.117)  IpH.TTH.TlTA 

34.  WHITE  (4,114) 

33.  DO  107  KH«1.N 

34.  btllTL  (6.119)  HK.SPRIKH).SIR(HK).SCTAIKK) 

37,  107  CONTINUE 

36.  104  CONTINUE 

39.  STOP 

C 

60.  109  FONHAT  I1H1.42X.31H  PUHOUC  UNIVrHSllT  COhlPCSSGR  PROOHAh  * PUCPP 

1*001/////////) 

61.  110  FOHHAT  (213) 

62.  Ill  FOHhAT  ( 2F  10 .0  ) 

63.  112  FOHhAT  (3F1U.0) 

64.  113  FOHHAT  113) 

63,  114  FOHhAT  I6F10.0) 

66.  113  FOPHM  (////1I.6H  SIAGL.l2.32l.  OAIA  IS  NOT  SUFFICIENT  - CALCULATIO 

IN  TLHHINATEU  FUH.F6.4.10H  HASS  FLOW) 

67.  116  FOI'HAI  1////1I.13X  HASS  Flow  I . F 1 0 . 4/2X  . 6HHPP  s .F10.4/1X  ,27m  INL 

ID  TOT  AL  TthhIHAIuHL  = .F  | 0 ,4/ 1 X . 24F,  INLtI  TOTAL  PRISLUPI  * .F10.4 
2/ 1 X . 26*.  INkL  I HOLfCULAK  WrlGMl  = .F  10.41 
66.  117  FOHhAT  (/IX. 24m  COhPHLSSUH  PMESSUHl  PATIO  = .F 1 0 , 4/1 X , 32H  CohPRCSS 

10H  IIHIIHATUKL  PATIO  « .F  10.4/1X.33H  COF.PPISSGP  AU 1 AHA  TIL  LFFIClEFt 
2CT  A . F 10,4* 

69.  116  FoPhAT  I//I0A.21*.  STAGE  PhLSsOPT  PaT10.10X.24H  STAGE  TEhPEHATuRC  6 

lATIO.lbI.27H  STAGE  AUlAl.aflC  (FFUUNCT) 

70.  119  FOHhAT  ( /2A.6I1S1A(,L  , 1 1 . Ox  .F  1 0 . 6 . 2 3A  . F 1 0 , 6 . 24  X iF  10 ,6) 

C 

Tt.  ENO 
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APPFNMX  C 


PURDUE  UNIVERSITY  COMPRESSOR  PERFORMANCE  PR0GRAM-PUCPP-002 

During  the  analytical  examination  of  the  N.G.T.E.  Reference  #109 
Compressor  it  was  desired  to  determine  the  effects  of  mass  flow  bleed 
and  injection  toward  negating  the  performance  change  caused  by  the 
liquid  phase  evaporation.  This  program  was  developed  to  specify  the 
amount  of  mass  flow  bleed  and/or  injection  needed  to  return  the  com- 
pressor to  the  point  which  give  the  same  pressure  ratio  as  its  original 
operating  point.  The  analysis  also  includes  the  possibility  of  heating 
or  cooling  the  injected  gas,  as  desired  by  the  user.  The  following 
four  sections  describe  the  necessary  input  data,  the  resulting  output 
data,  the  method  and  capabilities  of  the  program  and  the  program  listing. 

1.  Input  Data:  There  are  two  items  of  input:  (i)  the  corrected 
performance  data  at  the  operating  point,  and  the  necessary  inlet  condi- 
tions used  to  obtain  their  corrected  forms,  and  (ii)  factors  that  govern 
the  injection  process.  Both  integer  and  real  input  formats  are  used. 
Integer  numbers  are  punched  in  fields  of  two  locations  starting  with 
column  1.  No  decimal  points  may  be  used.  Real  numbers  are  punched 
in  fields  of  ten  locations  starting  with  column  1.  Decimal  points  must 
be  used  for  real  numbers  if  right-hand  justification  in  each  field  is  to 
be  ignored. 

A listing  of  the  input  cards  is  given  in  the  following  chart, 
where  one  line  represents  one  card. 
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Note  that  Figure  2.9.  of  Report  No.  M-WPAFB-77-1  which  is  referred 
to  here  is  included  at  the  end  of  this  Appendix  as  Figure  6.1. 


N 


WCl 

WCA 

RC1 

RCA 

RPM2 

Occurs 

W1 

T01 

P01 

FL0W1 

WI 

TO  I 

N 

Times 

CP 

DELTP 

DELTN 

TMAXP 

TMAXN 

XMAX 

YMAX 

Variable  definitions  - The  input  data  must  have  units  corresponding 
to  the  dimensions  given  below  (were  applicable). 

Length  feet 

Time  seconds 

Mass  pound-mass 

Force  pound-force 

Temperature  Degrees  Rankine 

The  one  exception  involves  the  units  of  RPM2,  which  are  per  minute. 

N Number  of  calculations  to  be  made. 

WC1  Corrected  mass  flow  of  point  1 (see  Figure  6.1) 

Corrected  mass  flow  of  point  a (see  Figure  6.1) 

Corrected  rotor  speed  of  point  1 (see  Figure  6.1) 

Corrected  rotor  speed  of  point  a (see  Figure  6.1) 


WCA 

RC1 

RCA 

RPM2 


Compressor  speed  at  point  2.  This  may  be  arbitrarily 
chosen.  Unitsiminute"1 . 


W1  Molecular  weight  of  conditions  at  point  1 (see  Figure 

6.1 ) Uni ts :pound-mass/pound-mole. 

T01  Total  temperature  of  conditions  at  point  1 (see  Figure 

6.1) Units:Deg.  Rankine. 

P01  Total  pressure  of  conditions  at  point  1 (see  Figure 

6.1 ) Uni ts :pound-force/feet2 . 

FL0W1  Mass  flow  at  point  1 (see  Figure  6.1) 

Uni ts: pound-mass/second 


I 

I 

i 


4 


1 
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WI  Molecular  weight  of  gas  to  be  injected. 

Uni ts :pound-mass/pound-mole. 

TOI  Total  temperature  of  gas  to  be  injected,  before 

heating  or  cooling  occurs.  UnitsrDegrees  Rankine. 

CP  Specific  heat  at  constant  pressure.  Assumed  to  be 
constant  throughout  compressor. 

Uni  ts :feet-pounds-force/pound-mass/Degrees  Rankine 

DELTP  Specifies  the  incremental  change  in  temperature  of 
heated  injected  air  that  is  desired  during  the 
analysis.  For  a given  mass  flow,  increments  of 
DELTP  will  be  added  to  the  injection  temperature  and 
analyzed  until  TMAXP  is  exceeded. 

Uni ts:Degrees  Rankine. 

TMAXP  Maximum  temperature  for  heated  injected  gas.  Note 
that  TMAXP/DELTP  must  be  an  integer. 

Units: Degrees  Rankine 

DELTN  Specifies  the  incremental  change  in  temperature  of 
the  cooled  injected  gas  that  is  desired  during  the 
analysis.  For  a given  mass  flow,  increments  of 
DELTN  will  be  subtracted  from  the  injection  temper- 
ature and  analyzed  until  TMAXN  is  exceeded. 

TMAXN  Minimum  temperature  for  cooled  injected  gas.  Note 
that  TMAXN/DELTN  must  be  an  integer. 

Units :Degrees  Rankine. 

XMAX  Maximum  injection  mass  flow  ratio.  This  eliminates 
physically  unrealizable  bleed  rates. 

YMAX  Maximum  bleed  mass  flow  ratio.  This  eliminates 
physically  unrealizable  bleed  rates. 

Note  that  TMAXP/DELTP  + TMAXN/DELTN  £ 20. 

To  insure  correct  non-dimensional  variables  being  input,  the  following 
equation  should  be  used. 


T 


where,  tfi  = 


N = 


WC  = 
RC  = 


( ) 


x 


mass  flow 

rotor  rotational  speed 
stage  inlet  total  temperature 
stage  inlet  total  pressure 
corrected  mass  flow 
corrected  RPM 

pertaining  to  the  operating  point  x (x  = a,  1) 


2.  Output  Data:  There  are  two  types  of  output  data:  (i)  the  injection 
and  bleed  mass  flows  for  the  given  operating  conditions,  and  (ii)  diag- 
nostic explanation  concerning  physically  impossible  developments. 

Providing  the  latter  does  not  occur,  the  program  will  give  the  follow- 
ing results: 

1.  Injection  mass  flow  ratio,  x 

2.  Bleed  mass  flow  ratio,  y. 

In  addition,  the  input  data  and  injection  conditions  corresponding  to 
these  input  data  will  be  printed.  If  values  of  x and/or  y are  negative 
or  greater  than  the  maximum  limits,  as  given  by  XMAX  and  YMAX  respectively, 
the  following  statements  will  be  printed: 

VALUES  OF  X OR  Y ARE  NEGATIVE 
VALUES  OF  X OR  Y EXCEED  LIMIT 
respectively. 

3.  General  Equations  and  Method:  Figure  6.1  illustrates  the  desired 
result  of  returning  a certain  operating  point  to  the  position  which 
give  the  same  pressure  ratio  as  its  original  position  by  changing  the 
corrected  parameters  associated  with  the  point.  The  method  of  altering 


i 


60 


the  parameters  is  the  bleeding  and  injection  at  particular  stage  inlet. 
If  point  (1)  (in  Fig.  6.1)  is  to  return  to  the  point  (2)  which  give 
the  same  pressure  ratio  as  its  original  point  (a),  it  is  necessary  that 
the  following  expressions  be  satisfied. 


(SP)a  = (SP)2  C-3 

(MP)a  = (MP)2  C-4 


where  the  subscript  a and  2 corresponds  to  the  point  (a)  and  (2)  in 
Figure  6.1,  respectively,  and  this  will  hold  throughout  this 
analysis  and  where 


SP  = 

MP  = 


mw  = 


Y = 

M = 


Speed  Parameter  = 


+ M2)1/2 


m/T 

Mass  Flow  Parameter  = -p- 

0 


+ 


Y + 1 


Total  Temperature 
Total  Pressure 
molecular  weight 
specific  heat  ratio 
Mach  number 


i 


i 
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That  is 


2 


C-10 


B = 


!Vi 


n 


l + 


Y,  - 1 


M2 

a 


1 + 


Y9  - 1 


— M2 


C-ll 


n 


Ya  + 1 
y - 1 y - 1 
(1  + -^—M2)  a 

y2  + 1 
Y?  " 1 Yo  “ 1 


C-12 


In  deriving  the  equations  ( C-7  ) and  ( C-8  )»  it  was  assumed  that  the 
pressure  is  unaffected  by  small  amounts  of  injection  and/or  bleed  and 
the  following  relationship  was  used. 


m2  = m1  (1  + x-y) 


C-l  3 


where 

_ minject 

A “ * 

ml 

a mbleed 
y ~ m. 
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The  molecular  weight  of  mixture  after  injection  and/or  bleed,  mw2. 


can  be  determined  as  follows. 


m,  m.  m, 

mw0  = mw,  + mw.  -r*-  - mw.  -r— 
2 1 m2  i 1^2  b m2 


where 


m, ,mw,  = mass  flow  and  molecular  weight  of  mixture  corresponding 
to  point  ( 1 ) 

iti?,mw?  = mass  flow  and  molecular  weight  of  mixture  corresponding 
to  point  (2) 

itu ,mwi  = mass  flow  and  molecular  weight  of  mixture  to  be  injected 
^b,mwb  = mass  ^ow  anc*  mo^ecu^ar  ^ight  of  mixture  to  bleed 


Assuming  mwb  = mw1 , Eq.  (A. 3. 14)  becomes 
(1  - y) .mw,  + x.mw. 

"“2  ■ --IT,-) C-,S 

In  the  same  manner,  the  total  temperature  corresponding  to  point  (2),  , 

is  determined  by 


(1  - y)  t01  + Z-X-Tpi 

1 + x-y 


where 


T0i  = total  temperature  of  mixture  to  be  injected 

Tq1  = total  temperature  of  mixture  corresponding  to  point  (1) 

Z = multiplicative  factor  relating  the  heated  or  cooled 
injection  total  temperature  to  Ti 

^heated^i  or  Tcooled^Ti^ 
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Assuming  that  the  specific  heat  at  constant  pressure  is  constant  through- 
out compressor,  the  specific  heat  ratio  corresponding  to  point  (2), 

Y2,  is  given  by 


y2  = 


1 - 


Ru 


C .mw0 
P 2 


C-l  7 


where  Cp  = specific  heat  at  constant  pressure 
Ru  = universal  gas  constant 

Substituting  equations  ( C-l 5 ),  ( C-l 6 ),  and  ( C-l 7 ) into  equation 
( c-7  ) and  ( C-8  ).  and  solving  for  x and  y,  one  can  obtain 

r (mw.  + C.mw, ) 

_ _ t-  1 * r ip 


y = 1 - C x 


C-19 


where 


C = 


zTq.A  - mw.  + Ru/Cp 

— T o i A + mw1  - Ru/Cp 


A = 


E - !Vi/- H. . m 

m,  V Ya.mwa 


C-20 


C-21 

C-22 
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(1  + C)] 


C-23 


D = (zTi  + CT] ) [mw.  + C.mw1  - g 


Ru 


F = /D 


C-24 


,1 


Equations  ( C-18  ) through  ( C-24  ) are  utilized  in  the  program  to 
directly  determine  values  of  x and  y.  It  is  also  of  interest  to  deter- 
mine the  values  of  x and  y for  the  cases  of  injection  only  and  bleed 
only,  respectively.  At  present,  only  the  case  of  bleed  only  has  been 
incorporated  into  the  program.  For  this  condition,  the  following 
equations  must  be  satisfied. 


C-25 


C-26 


In  the  present  case,  the  value  of  £ and  n.  which  are  given  by  Eqs.  (C-ll 
and  C-12  ),  are  considered  to  be  nearly  one.  Thus,  assuming  that 
€ = n = 1 , the  proqram  was  developed  to  specify  the  amount  of  mass  flow 
bleed  and/or  injection  needed  to  return  the  compressor  to  the  point 
which  give  the  same  pressure  ratio  as  its  original  operating  point. 

4.  The  Program  Listing: 
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1.  dimension  zizii.  mm.  ti2u 

2.  INTEGER  TINCP 

3.  INTEGER  TINCN 

k.  I'ilSGS.O 

5.  .HlTE  (4.1081 

4.  108  FORMAT  IlHl.kOX.SlH  PUROUE  UNIVERSITY  COMPRESSOR  PROGRAM  - PUCPP 

1-002///////) 

T.  R) »0  (S. 1091  N 

8.  10*  Format  ( 12) 

».  OO  107  1=1. N 

10.  RlAUIS.110IWCl.WCA.RCl.RCA. RPM2 

11.  110  FORMAT  (SF10.0I 

12.  Rf  AO  (5.1111  Wl.Tol.POl.FtOWl.Wl.TOl 

IS.  Ill  FORMAT  (4F10.0I 

tk.  READ  (5.112)  CP, DrL  TP .OECTN, TMAXP .TMAXN.XMAX.YMAX 

IS.  112  FORMAT  (7F10.0) 

14.  WRITE  14.  US)  WCl  .WCa.RCi.RCa.RPM2.W1,T01.P01  .FLOW1.WI  .TOI  .CP. 

10FLTP. DClTN, TMAXP, TmAXN.XrAX.VMAX 

17.  115  F0RMAT(////1X. 'INPUT  OATA'/IX. »wClc' .FlO.k.SX. *WCAs' .FlO.k.SX, 

2'RC1=' .F10.R.5X, 'RCA= 'FlO.k.SX, 'HPM2=' .F 10 . k , SX/l X . 'Wl»' .FlO.k, 
35X. 'T01=' .FlO.k.SX. >P01=' .FlO.k.SX, 'FL0W1 *• ,F* .k , 5X/J X . ' W 1= • ,F9,6, 
RSX. *TOI=' .FlO.k.SX. 'CPs' .FlO.k.SX, 'DEL TP=' . 

SFlO.k.SX. 'DELTM*' .FlO.k. Sx/1X.'TMAXPe', FlO.k. SX.'TMAXNs' 'FlO.k.SX, 
4'XMAX=' .FlO.k.SX.'YMAXs', FlO.k.SX///) 

18.  UNCG=1 .0 

1*.  2(ll»l.O 

20.  IF  (AI)S)0£LTP),LT..0001)  OELTPsI.O 

21.  IF  (ADS(0CETNI.lT..0001)  OElTN>1.0 

22.  TINCPSTMAXP/OELTP.1.0 

23.  TINCNsTMAXN/UELTN 

2k.  As(RCA/RHM2)**2k(28.96k/l,koi 

2S,  E=(WCA»P01/FL0W1 )/S0RT 1 1 . kO ) /SORT (28,94k) 

24.  KRSTINCN.TINCP 

27.  UELTsOEETP 

28.  00  102  K=1,KK 

29.  IF  (R.EO.l)  GO  TO  101 

30.  2(K)s(T0I*UNEG«(lCLT)/T0I 

81.  DELT=DELT»UELTP 

32.  101  Cs|Z(H)«A«T0I-WI.R/CP)/(Wl-A*T01-R/CP) 

33.  O=(Z(K)*TOI»C*TOl)*(Wl*C*wl-R/CP*(1.0.C)) 

3k.  FsOss.S 

33.  X(R)sE/Fk(WI«C>wl)/(l.O«Ct 

34.  T(K)zl,0-C*X(K) 

37.  IF  (K.NE. TINCP)  GO  TO  102 

38.  bNEGs-1.0 

39.  OELTPSUELTN 

kO.  OELT  sOELTN 

kl.  102  CONTINUE 

k2.  RPH22=RCA/|Wl/T01-R/|CP»T0l> > *>0 . SkSORT I 28 .*4k > /SORT  1 1 .k ) 

k3.  YY  = 1.0-E/(T01/WI-R/CPM01/(W1)**2I**,S 

kk,  T022=T01 

kS.  FLOW22sFLOWl*(1.0-YY) 

4*.  P02=P0l 

kT.  HC22=RPM22/S8R)(T022) 

k8.  WC22sFL0w22/P02«S0RT(T022) 

k*.  IF  ((ABS(RPM2-RPM22) /RPM2 ) • LE • • 001 1 GO  TO  104 

SO.  00  10S  Lsl.RR 

91.  T 0=Z 1 L ) *T0 1 

32.  IF  (X(l I.GC.O.O.ANO.Y(LI.GE.O.O)  GO  TO  103 

53.  WRITE  (4,115)  TO 

SR,  115  F0RMAT(//1X. 'TUAOn*' .FlO.k, 18X.30H  VALUES  OF  X OR  Y ARE  NEGATIVE) 
99.  GO  TO  105 

54.  103  IF  ( X ( L I .LE . XMAX . A NO . Y ( L ) .LE • YMAX ) GO  TO  10R 

37.  WRITE  (4,114)  TO 

98.  114  F0HMaT(//1X«  ' TUAOOk ' , FlO.k, 18X,  'VALUES  OF  X OR  Y EXCEED  UNIT  •) 

3*.  GO  TO  109 

40.  10k  T02=l ll.O-YIL) ) «T01 ♦ T0» X ( L ) I / ( 1 . 0«x ( L » -Y ( L I I 

41.  FLOW2sFLOWl*(1.0kX(L)-YIU I 

42.  P02=P0l 

43.  RC2»RPM2/SURTITn2) 

4k.  WC2=FL0W2/P02«S0RT( T02I 

49.  WRITE (4, 117)  TO.X(L).YIL) .T02.FL0W2.P02.RC2.WC2 

44.  117  FORMAT ( //IX, 'TOAOOs. ,F9.k.3X. 'X«' ,F7,k,3X, 'Yk'.FT.k.SX, 'T02»« ,F9.k 

1 . 3X , • FI  Ow2= • ,F9,k,3X.'P02*'.F*.k.3X.'RC2s' .Fk . R . SX . <WC2« ' .F*,k 1 

47.  103  CONTINUE 

48.  IF1YY. LC.YRAX.AN0.YY.6C, 0,01  GO  TO  121 

4*.  WR ITE ( 4 , 120  t 

TO.  120  FORMAT(//IX, 'THIS  CANNOT  BC  ACC0NPL1SHE0  WITH  BLEEO  ONLY') 

71.  GO  TO  107 

72.  121  WR!TC(4.U*>  RPN27.YY.RC22.WC22 

73.  118  FORMATI/ZIX.'FOR  NLFEO  ONiY'.3x.'RPM*'.F10.k,5X.'Y»',F7.k,9x. 

2 • RC2* • .FlO.k.SX, ,wC2*' .Flo. 9) 

7k.  GO  TO  107 

73.  104  WRIT£(4«119)  RPM22.YY.RC22.WC22 

74.  11*  FORMAT (//IX.'THiS  MAY  BC  ACCOMPLISHED  WITH  BlCEO  ONL Y« . 9X . 'RPm« • . 

SFlO.k.SX. 'Ys* ,F7.k.SX,'RC2*' ,F 10 . * . SX , • WC2» • , F 10 . 5 ) 

77.  107  CONTINUE 

78.  STOP 


7* 


C 

C 


End 
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Appendix  D 

EQUATIONS  OF  MOTION  FOR  TWO  PHASE  FLOW 


I 


Nomenclature 

Reference  should  be  made  to  Report  No.  AFAPL-TR-76-77  for  the  nomen- 
clature except  as  noted  below. 

Ej  internal  energy  per  unit  mass  of  the  class  i droplets  ( i =1 ,2) 

E energy  required  in  connection  with  forces  stem  from  momentum 

v change  because  of  interaction  between  the  two  classes  of  droplets 

E energy  required  for  the  change  of  droplet  surface  area  during 
formation  of  new  droplets 

m mass  that  would  be  created  due  to  condensation  per  unit  volume 
of  mixture 

rh  mass  that  would  be  lost  due  to  evaporation  per  unit  volume  of 
mixture 

m.j  average  mass  of  a droplet  of  class  I 

N.  total  number  of  droplets  generated  per  unit  volume  per  unit  time 

due  to  break  up 

ag  bulk  volume  concentration  of  drizzle  in  the  mixture 

ij>.  forces  exerted  on  the  phase  of  class  i droplets  by  the  other  class 
of  droplets  ( i=l ,2) 


I 
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I.  Three-Dimensional  Conservation  Equations  in  Intrinsic  Coordinates 


(a)  Mass  conservation  Equations. 

Gas  Phase: 

+g])+£( 

+d)  +£(U>  f(i  - +fi])  + £(  rU4(  1 ~ v-*)/? 

+/>])]  = ('*«.  ■'*«.)  + *■<:*)  ( D-l  ) 

Liquid  Phase: 

1 . Droplet  - Type  1: 


2.  Droplet  - Type  2: 

t[w(r£K»J+£(£^)+&(rfi»/r,)+£;(r/2Ur*) 

+ £(/2Me)-,-£r(l'/iU-»)] 

= (orica  - ) -1-Mi  “ N|  ^1  ( D-3  ) 
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(b)  Momentum  Conservation  Equations. 


Gas  Phase: 


m-component 

[( ' -*-*)/?  +/?]  [( L4. ) (#“  - #) 


il 


(D-4) 


9 -component 


[('-«?  - *)/?  +/?]  F t(  ^ U-)(  gp* * f^) 

►rK  + u)(#+#)J 


(D-5) 


n-component 

= " ^ + + ^ +'Fi»  (D-6) 
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< 

t 

i 


4 


Liquid  Phase: 


1 . Droplet  - Type  1 : 

m-component 

-fL-*-  t,  + N,  ) 

- M,  ( lA^m,  + Un  ) (D-7) 


9 -component 

4 M UJ(-§^  +£r ) + r(K,-U, ) • 

+ N/|  ™,  (\A/Fetz  + Ue)-  M,(V\/i,o,->-Ue)  (D-8) 


n-component 


+ M (Wr2  + y„)-M,(Wr,  + Ur>)  (D-9) 
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2.  Droplet  - Type  2: 


9 Wpmj  9 Lh*  ^ "] -7—  -1 — 1 — , 

+ M,(l^m|+U-n,)--Nl  ^l(^+Uw) 


(D-10) 


6-component 

ff[H  ) + r ■ 

+ u»)-w  m(W£„1-U()  (D-11) 

n-component 

/? [( Mf,  +U.) f-#2  + #)  - #M£.+ 14 )( ( li> 


3 Mpru  _ 9 Un 

**  Be  Be 


(W^  + UJ-N/!  =wr,  (W,.,+U„)  (0-12) 


(c)  Energy  Conservation  Equations. 


Gas  Phase: 

d§*+lr)]=  $ii‘fI  + ^ + ^-*-w„)^'+7-(^ 


+ W*) 


i£  _ JIta.  <9^ 

90  <?777  3777 


-dfrf  + d/jg 
90  + a/9 


) 


- -g.far’.  _ 

an 


( D- 13) 
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Liquid  Phase: 


1 . Droplet  - Type  1: 

v\{ %.,+U.)-^1  + = ij  + ¥f, 

_ dh*»>  _ dffr-mi L(  Jl$hi  _ 

9th  ?*i  n 9^  / «??» 

+ # - M,  Ej,  +M  ^EXj  -+ES,  - E„  - E„ 

(D-14) 


2.  Droplet  - Type  2: 

dip**  , J_f  dftpft?  _ dfriez  A _ 99^,2 

?7n  ~ dm  ' I'l  3^  <90  / 5>7i 

+ 4^  + w,Ea-E,,^-F» 

(D-l 5) 


II . Axisymmetric  Conservation  Equations  in  Intrinsic  Coordinates 


(a)  Mass  Conservation  Equations. 

Gas  Phase: 

7 f^T ( t [( ' - V" - <4 )f ■*•«] W<-)  + ■ t [( '■ ~ -■ *)%■ +/?]lr] 
=(™<,~n,)+-<  •*%)  (D-i6) 
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Liquid  Phase: 

1 . Droplet  - Type  1: 


(D-17) 


2.  Droplet  - Type  2: 


(D-18) 


I 


(b)  Momentum  Conservation  Equations. 


Gas  Phase: 
m-component 


9 -component 


(D-20) 


n-component 


[( i - *-*>**]  [ -H  * ■ - ru>jV£ J 

= “i£  + ^,+^  + 7i'„ 


(D-21 ) 
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Liquid  Phase: 


m-component 


1 . Droplet  - Type  1 


=Tj  +r  -Ffc 


+%,  + n,™,(  Tft'+iA*)  - n (w,~.  +u«) 


(D-22) 


0-component 


t-’H -1^+^  -7-,  - r„ 


+ M ■U0)-Ml(]A/pei-i-Uff) 


(D-23) 


n-component 


+ t + M (Wp„+  U„)-  M,  W.) 


(D-24) 


2.  Droplet  - Type  2 : 


m-component 


n ("w  _i 

'2  L ^ r 


+ t,  + M‘  (n£-.  i-  u.)  - w,  f w^,2-h  u„) 


(D-25) 


0-component 


a -^t^r  ( r i%. ) - £ 

4*  4^  4-  Mi  ( lA/p«»i  ■♦•Me)  - IS/,  ( ^Vp«  + Me) 

(D-26) 
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s[w. 


n-component 
3 M/p -n  a 

j»tni  97n 


+ ^ -+-  M,  ( 1 Vr,  ■+•  U*  ) - /V,  Tn,  ( kVp^i  -+-  U-n) 

(D-27) 


(c)  Energy  Conservation  Equations. 
Gas  Phase: 


[(•-*-*)£  +/»]  M4.  (#£  + |£)  = %+T? 


4.  ^ ^ ^ _ 2#Tt, 

8'm  <?W1  <?>77l  377  3n 


( D-28) 


Liquid  Phase: 


1 . Droplet  - Type  1: 


p W =S  $ - $ 4-  b - — legt  4-  J&2<  - .glidi 

'•  dm  *pi  Xj,  ^ ^ dm  dn 


+ j*h  - m,eIi+w,w,eIi+fj,-e„-f, 


n 

(D-29) 


m 


du,. 

f*"2  2-771 


2.  Droplet  - Type  2: 

= $ - s ~t-~F  - Jhrt  4.  ~ l£ei 

-‘y2  *2*  ip,  5>7n  3771  3 -n 


9f 


+ + M'  El,  - ^-5.  *5*  - Efc 


( D-30) 
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Ill . Three-Dimensional  Conservation  Equations  in  Cylindrical  Coordinates 

(a)  Mass  Conservation  Equations. 

Gas  Phase: 

- (-We,-  •wic,)  + ('«'ei-  -”<cj  <D_3n 


Liquid  Phase 

1 . Droplet  - Type  1 : 


- fac,  - 'We,)  — M,  + M (d-32) 


2.  Droplet  - Type  2: 


= (™c*  ~ 771*  ) -h  M,  “ A/j  sm,  ( D-33) 


(b)  Momentum  Conservation  Equations. 
Gas  Phase: 


r-component 


( D-34) 
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79 
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(c)  Energy  Conservation  Equations. 


Gas  Phase:  , 

[0-<v -<*►)/?  ^][^r(TF 


dht>  \ . f 

zir  )^~{  r 


-h  CO 


+ ^ ^ f "T^  + “O'SF  ’¥'w9t~Ti- 

3%r  _ dfr>- I r 

?r  ~ r l 20  30  J ~ dz  ~ bz 

( D-43) 


Liquid  Phase: 


1 . Droplet  - Type  1 : 


* [%.  If*  ' + (■ ^ + *>)  & + %,  If*  ] = ih  ■ - f*  + 1 

— dftrw  J_  f dfo»i  3 'l  3^pti  , 

a/-  ~ r\30  3&  ) a?  ar 

+ -Elk  - m.F*  + M ^ E*  + Esi  - E„  - Ew 


/?[% 


( D-44) 

2.  Droplet  - Type  2: 

i!te  + (- 

ar  v 

14^*2  > 

"F” 

l-w^+Mfc, 

3?  J ^ iTa 

<^#>2 

1 / 3^2  _ 

r 7* 

ar 

rl  a* 

33  / 5>Z 

fin 

?r 

Ml  E^'+’Hsa 

-E»- 

E>2 

( D-45) 
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(a)  Mass  Conservation  Equations. 


Gas  Phase: 

= (^V1-^e,)+-C  '*Vr2-  ™Cj.)  ( D-46 ) 


Liquid  Phase: 


1 . Droplet  - Type  1 : 


r dr  ( Yf + 31  0?  ^n)  ~ fA:,-  ™e,)-  M,+A/,  m, 

(D-47) 

2.  Droplet  - Type  2: 

7 w(  Yfr\ /Vrz)+-k(£  Wfa)  ss('*ea“*hek)  + MrMi 


( D-48) 


(b)  Momentum  conservation  Equations. 

Gas  Phase: 
r-component 


?/* 


(D-49) 


6-component 


[( i-  *-<*  )yf  +^]  f (-1^  + «>)  -r  %.  (-&*+ U>) 


(D-50) 
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(D-51) 


Liquid  Phase: 


4 


r-component 


1 . Droplet  - Type  1: 


W, 


BlVP 


m _ 


( + ru>)z 


+ woi 


dH 


m 


J 


frt  y~  ' 'P*1 

fc  + N,  Wrn  ~ nWrn 


(D-52) 


6 -component 

= +£P„-  ■&  + 


-Mi  ( Wt.,+ 


(D-53) 


z-component 


ar 


an£ 


Wi 


az 


(D-54) 
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2.  Droplet  - Type  2: 


r-component 


4% 


frz  a r 


( r^)  / 

air  ~ +"  gg 


’Tip ■*■  Mi  M^r.  — N/,  ^ 


(D-55) 


0-component 


/?K«  (-§7^ + “>)+%*  (-^fe+ «.)+■  ^ 

= + ¥p«  " !}«  ■+■  fn  + Mi  ( ruJ) 

-Ml  W,  + lruj)  (o-56) 


z-component 


4 %*  T^5  + ^-#2J=  Tfp„+ V - ^ + fc 

+ Mi  l^Vp^,  ~ N/j  O^rT,  W£Zz  (D-57) 


(c)  Energy  Conservation  Equations. 
Gas  Phase: 


[('-*-*)/?  +/j]  [ (# +#)+nfc  (#  +#’)] 


= $3  + §1+H+Wr-0-+%4r--0c 


^ fijr  9 9 

~ ~ dz  dz 


(D-58) 


5 


Liquid  Phase 


1 . Droplet  - Type  1: 


i A/  •+■  M/  4^-1=  $ - $ 4- 

'*  L KVpn  ar  wpgl  a z J xpi  xi, T rh  ^ ar 


9 <b?3\  , 3 


3 Z 


+ - M,Et,  + N,  E^+Ej, 


- ES1  - E 


fl 


(D-59) 


2.  Droplet  - Type  2: 


/?[>,„ Ih-  §*  + ?„  - 


30 


m 


dir 


3 %.\r3  3 <5^X« 


, grira  o ^ t-  _ » / -r- 

+ ar  az  + dz  M 


+ Ett~  Ej,  -E 


't'z 


(D-60) 
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Appendix  E 

MODIFICATIONS  TO  UD0-300  TO  ACCOUNT 
FOR  VARIABLE  MOLECULAR  WEIGHT,  SPECIFIC 
HEATS  AND  ACOUSTIC  VELOCITY 

The  program  UD0-300  utilizes  a series  of  FUNCTION  subprograms  in 
determining  all  thermodynamic  and  flow  properties  involving  molecular 
weight,  specific  heat  at  constant  pressure  and  acoustic  velocity.  Each 
FUNCTION  subprogram  corresponds  to  a specific  gas-phase  property.  These 
subprograms  may  be  separated  into  three  distinct  categories,  each  of  which 
must  be  handled  differently. 

1.  FUNCTION  Subprogram  UDG1 : used  to  input  certain  gas-phase  properties 
to  UD0-300.  In  the  original  program,  UDGl  input  the  constant  value  of 
specific  heat  at  constant  pressure,  and  molecular  weight  in  the  form  of  a 
gas  constant.  Since  both  these  properties  are  variable  for  two-phase  flow 
considerations,  UDGl  will  now  input  the  compressor  inlet  values  of  specific 
heat  at  constant  pressure  and  molecular  weight.  These  will  serve  as  refer- 
ence values.  Also,  certain  properties  pertaining  to  the  particulate  phase 
will  be  input  by  this  subprogram  for  use  in  the  calculation  of  the  gas-phase 
acoustic  velocity. 

2.  FUNCTION  Subprograms  UDG3-UDG5  and  UDG7-  UDG9:  output  certain  gas  pro- 
perties using  input  of  enthalpy.  Each  of  these  subprograms  determines  a 
result  from  an  equation  of  the  form: 

x = 4>(c  ) f(mw  ) 6(H)  g(?) 
y 
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where,  x 
4>(c  ) 

Hg 

f (mw  ) 
6(H) 


desired  flow  property 
dependence  of  x on  cp 

dependence  of  x on  mw 


9 

g 


dependence  of  x on  enthalpy,  H 


g(?)  = some  arbitrary  function  unique  to  each  subprogram  providing 

equality  (may  involve  functions  of  additional  flow  properties). 


Using  the  perfect  gas  assumption,  that  the  gas-phase  specific  heat  and 
recalling  and  molecular  weight  are  given  as  functions  of  the  gas-phase 
mixture  temperature  only,  then 


c = H/T  = f(H/c  ) = f,(H)  (E-l) 

Hg  y Hg 

and 

mwg  = mwg(Tg)  = mwg(cp  ) = f2(H).  (E-2) 

An  iterative  technique  is  needed  to  determine  c and  mw  due  to  the  nature 

Pg  g 

of  equations  (E-l  ) and  (E-2  ), given  only  the  enthalpy.  An  addition  to  these 

i 

subprograms  that  is  needed  to  account  for  variable  specific  heat  and  mole- 
cular weight,  therefore  consists  in  the  iterative  technique. 

In  the  case  of  subprogram  UDG9,  g(?)  involves  the  acoustic  velocity. 

An  extra  section  is  added  to  this  subprogram  to  account  for  droplet  affects 
• on  the  gas-phase  speed  of  sound. 

3.  FUNCTION  Subprograms  UDG2  and  UDG6;  output  enthalpy  using  input  of  cer- 
tain gas  properties.  In  general,  an  iteration  scheme  is  employed.  Enthalpy 
is  determined  using  the  reference  value  of  c and  mw  , input  by  U0G1 . 

Pg  9 

Using  the  equations  (E-l  ) and  ( E-2  ) of  the  iteration  technique  involved 

in  the  previous  FUNCTION  subprogram  discussion,  the  values  of  c and  mw 

.9 

corresponding  to  the  just  calculated  enthalpy  can  be  compared  with  the 


* 


i 


< 

j 
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the  reference  values.  A value  of  enthalpy  is  then  output  when  the  iter- 
ation scheme  converges  on  the  gas-phase  specific  heat  and  molecular  weight. 

ANALYSIS: 

(a)  Variable  molecular  weight  and  specific  heats 

- Subroutine  ITER:  This  subroutine  is  called  by  FUNCTION  sub- 
programs UDG2  - UDG5  and  UDG7  - UDG9.  A value  of  enthalpy  is 
input  to  the  subroutine  and  corresponding  gas-phase  specific 
heat  at  constant  pressure  and  molecular  weight  are  determined 

by  an  iterative  approach.  These  values  of  c and  mw  are 

P9  9 

returned  to  the  calling  FUNCTION  subprogram  by  the  labeled 

common  "GAS".  It  is  assumed  that  the  functions  of  c „ . 

p and 
9 

mWg  given  by  equations  ( E-l  ) and  ( E-2  ) are  well  behaved, and 
convergence  always  occurs. 

- Subroutine  SET:  The  purpose  of  SET  is  to  retain  the  reference 

values  or  c and  mw  , input  by  UDG1 , that  correspond  to  the 

Pg  g 

compressor  inlet  conditions.  These  reference  values  areused 
in  the  first  iteration  performed  by  ITER.  The  final  iteration 
fixes  different  values  for  c and  mw  corresponding  to  the 

pg  9 

input  enthalpy  and  transmits  these  back  to  the  calling  FUNCTION 
subprogram  by  the  labelled  common  "GAS". 

- Subroutine  RETURN:  This  subroutine  accepts  the  reference  values 
of  gas-phase  specific  heat  and  molecular  weight  stored  by  SET 
and  re-establishes  these  values  in  the  labeled  common  "GAS". 
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Thus,  whenever  ITER  is  called,  the  following  procedure  occurs: 

(i)  Reference  c and  mw  values  are  retained  by  SET 

p q J 

a y 

(ii)  ITER  determines  c and  mw  corresponding  to  input 

Pg  9 

enthalpy  using  reference  c and  mw  in  first  iteration 

pg  ^ 

(iii)  ITER  returns  values  of  c and  mw  found  in  (ii)  to 

Pg  9 

calling  FUNCTION  subprogram  by  way  of  "GAS" 

(iv)  Subroutine  RETURN  resets  the  values  of  cp  and  mw^  in 
common  "GAS"  to  reference  values  stored  by  SET. 


(b)  Acoustic  Velocity 

Equation  (0-1)  and  the  liquid-phase  properties  pp  and  cp 
(introduced  to  the  subprogram  by  the  labeled  common  "LIQ",  and 
now  input  to  the  program  UD0-300  by  the  subprogram  UDG1 ) are 
added  to  the  subprogram  UDG9. 


Modifications  to  UD0-300 


FUNCTION  UDG1  (LNLT) 
common/GAS/G ,EJ ,R ,CP ,R0JCP 
common/l I Q/RHOP , AD 
LOG  1 = 1 
LOG  2 = 2 

READ  (LOG  1,  100)  CP,R,G,EJ 
READ  (LOG  1,  101)  RH0P,AP 

100  FORMAT  (4F12.0) 

101  FORMAT  (2F12.0) 

If  (CP , EQ .0.0)  CP  = 0.24 
If  (R.EQ.0.0)  R = 53.32 
If  (G.EQ.O.O)  G = 32.174 
If  (EJ ,EQ,0.0)  EJ  = 778.16 
WRITE  (LOG  2,110)  CP,R,G,EJ,RHOP,AP 
110  FORMAT  (2x,/,10x,37H  INLET  SPECIFIC  HEAT  AT  CONST 

1 PRESSURE, 2x,lH  = ,F8. 5,/ ,10x ,18H  INLET  GAS  CONSTANT, 21x, 

2 1H  =,F8.4,/,10x,22H  GRAVITATIONAL  CONSTANT,  17x,lH  =,F8.4, 

3 / ,10x ,17  H JOULES  EQUIVALENT, 22x,lH  =,  F8.3,/ ,10x,20H  LIQUID 

4 PHASE  DENSITY, 19x,lH  = ,F8.4,/ ,10x ,24H  LIQUID  ACOUSTIC  VELOCITY, 

5 1 5x , 1 H = ,F8.2) 

LNCT  = LNCT  + 5 
ROJCP  = R/(EJ  x CP) 

RETURN 

END 


FUNCTION  UDG2(S,P) 
common/GAS/G, EJ ,R, CP, ROJCP 
CALL  SET(RO,ROJCO,CPO) 

TOLRR  = .001 
I = 1 

1 H = CP*EXP(S/CP  + ROJCP*ALOGCP) 

CALL  ITER(H) 

If ( I . NE . 1 )G0  TO  2 
HO  = H 
1 = 1 + 1 
GO  TO  1 

2 If ( (ABS(H  - HO)/H).LE.TOLRR)GO  TO  3 
GO  TO  1 

3 UDG2  = H 

CALL  RETURN (RO,ROJCO,CPO) 

RETURN 

END 


FUNCTION  UDG3(P,H) 
common/GAS/G, EJ  ,R, CP, ROJCP 
CALL  SET(RO,ROJCO,CPO) 

CALL  ITER(H) 

UDG3  = CP*AL OG(H/CP)  - R/EJ*ALOG(P) 
CALL  RETURN (RO,ROJCO,CPO) 

RETURN 

END 


FUNCTION  UDG4(H,S) 
common/GAS/G ,EJ ,R ,CP ,ROJCP 
CALL  SET(RO,ROJCO,CPO) 

CALL  ITER(H) 

UDG4  = EXP(ALOG(H/CP)/ROJCP  - EJ/R*S) 
CALL  RETURN (RO,ROJCO,CPO) 

RETURN 

END 


FUNCTION  UDG5(H,S) 
common/GAS/G, EJ ,R,CP,ROJCP 
CALL  SET(RO,ROJCO,CPO) 

CALL  ITER(H) 

UDG5  = UDG4(H,S)/(R*H)*CP 
CALL  RETURN (R0,R0JC0,CP0) 
RETURN 
END 
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★ Install 


FUNCTION  UDG6(P,T) 
CP  = f (T) 

UDG6  = CP*T 

RETURN 

END 


FUNCTION  UDG7(H,S) 
common/GAS/G,EJ,R, CP, ROJCP 
CALL  SET(RO.ROJCO.CPO) 

CALL  ITER(H) 

UDG7  - H/CP 

CALL  RETURN ( RO , RO JCO , CPO ) 

RETURN 

END 


FUNCTION  UDG8(H,S) 
common/GAS/G , E J , R , CP ,ROJCP 
CALL  SET(RO,ROJCO,CPO) 

CALL  ITER(H) 

UDG8  = 1.0/ (1.0  - ROJCP) 
CALL  RETURN (RO.ROJCO, CPO) 
RETURN 
END 


I' 


FUNCTION  UDG9(H,S,V2) 
common/GAS/G,EJ ,R, CP, ROJCP 
common/LIQ/RHOP,AP 
CALL  SET(RO.ROJCO.CPO) 

CALL  ITER(H) 

GAMMA  = 1.0/ (1.0  - ROJCP) 

AG  = CP/(GAMMA*G*R*H) 

PG  = EXP (ALOG( H/CP) /ROJCP  - EJ/R*S) 

TG  = H/CP 
RHOG  = PG/(R*TG) 

SIGV  = f (H,S)  [-»  av  = av(H,S)] 

AMI  = (1.0  - SIGV)*RHOG  + SIGV*RHOP 

AM2  = (1.0  - SIGV)/(RH0G*AG**2.)  + SIGV/(RH0P*AP**2.) 

AM  = (AM1*AM2)**-. 5 

UDG9  = V2/AM 

CALL  RETURN(RO,ROJCO,CPO) 

RETURN 

END 
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SUBROUTINE  SET(RO,ROJCO,CPO) 
common/GAS/G , EJ  ,R ,CP , ROJCP 
RO  = R 

ROJCO  = ROJCP 
CPO  = CP 
RETURN 
END 


SUBROUTINE  RETURN(RO,ROJCO,CPO) 
common/GAS/U,EJ,R, CP, ROJCP 
R = RO 

ROJCP  = ROJCO 
CP  = CPO 
RETURN 
END 


SUBROUTINE  ITER(H) 
common/GAS/G, EJ,R, CP, ROJCP 
TOLRR  = .001 
11  = 1 

2 TGI  = H/CP 

3 CPI  = CPI (TGI) (Install) 

If (ABS ( (CP  -CP1)/CP1).LE.T0LRR)G0  TO  5 

If ( I .GT. 1 )G0  TO  4 

DELT2  = CPI  - CP 

CP  = CP/1.1 

TG2  = TGI 

I = 2 

GO  TO  2 

4 DELT1  = CPI  - CP 
DELTT  = TGI  - TG2 
TG2  = TGI 

TGI  = TGI  - DELTT/(DELT1  - DELT2)*DELT1 
DELT2  = DELT1 
CP  = H/TG1 
GO  TO  3 

5 R1  = R1 (TGI) (Install) 

TG  = H/CP1 

R = R(TG ) ( Ins tal 1 ) 

I f ( ABS ( ( R1  - R)/R).LE.TOLRR)GO  TO  6 
TOLRR  = T0LRR/2.0 
GO  TO  1 

6 ROJCP  = R/(EJ*CP1 ) 

RETURN 

END 
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